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FOREWORD 


This  report  was  prepared  by  the  Research  Division  of 
American  Optical  Company,  Southbridge,  Massachusetts,  under 
contract  Nonr  4875(00),  Optical  Inhomogeneities  in  Pumped  Lasers, 
for  the  Office  of  Naval  Research,  Washington,  D.  C.,  as  part  of 
Project  Defender . 

Principal  contributors  to  this  report  were  C.  J.  Koester, 
Introduction,  Review  of  Previous  Work  and  overall  compilation  of 
material;  E.  Snitzer,  Thermally  Stable  Cavities;  L.  W.  Smith, 
Calculation  of  Ttemperature  Profile;  S.  M.  Bergmann,  Effect  of  the 
Laser  Field  on  the  Index  of  Refraction  of  the  Laser  Rod;  and 
D.W.  Cuff,  End  Region  Stress  Analysis. 

The  report  is  published  in  two  volumes,  the  first 
unclassified  and  the  second,  a  classified  supplement  entitled 
Optical  Inhomogeneities  in  Pumped  Lasers(u),  Volume  II. 


ABSTRACT 


This  report  delineates  efforts  expended  during  the  one 
year  period  of  contract  Nonr  4875(00),  entitled  Optical  Inhomo- 
geneiti^s  in  Pumped  Lasers,  which  required  an  analysis  in  depth 
of  the.  ,.:dl  distortion  of  wavefronts  in  large  laser  reds  and  discs. 
TWo  approaches  for  elimination  of  this  distortion  were  considered: 
(1)  Establishment  of  a  uniform  temperature  throughout  the  laser 
during  pumping  and  (2)  Development  of  athermal  glass  parameters. 
Elements  of  both  approaches  wore  found  necessary  to  solve  the 
problem.  In  separate  sections  of  the  report  detailed  treatments 
are  given  to  Thermally  Stable  Cavities,  Temperature  Profile 
Calculations,  the  Effect  of  the  Laser  Field  on  Index  of 
Refraction  of  the  Laser  Rod,  and  finally  an  Analysis  of  End 
Region  Stress. 

A  classified  supplement  published  separately  as  Volume  II 
of  this  report  completes  the  research  findings  under  this  project. 
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1.  INTRODUCTION  and  SUMMARY 


The  effort  in  this  program  has  been  a  theoretical 
analysis  of  the  thernal  distortion  problem  in  glass  laser  rods 
and  discs.  The  goal  has  been  not  merely  analysis,  however,  but 
rather  a  more  thorough  understanding  of  the  causes  of  wavefront 
distortion,  leading  to  more  effective  athermalization . 

There  are  two  basic  approaches  to  eliminating  wavefront 
distortion  due  to  thermal  effects:  (1)  Attempt  to  establish  a 
uniform  temperature  throughout  the  laser  during  pumping . 

(2)  Develop  a  glass  in  which  the  change  in  length  due  to 
thermal  expansion  is  balanced  by  the  (negative)  change  in  index 
with  temperature .  The  nature  of  the  cavity  must  be  considered  in 
arriving  at  the  glass  specifications.  Also  the  effects  of  stress- 
birefringence  must  be  removed  or  reduced  to  an  acceptable  level. 

For  a  number  of  reasons  it  is  desirable  to  use  both 
approaches.  The  first  method  by  itself  will  probably  not  be 
sufficient  since  it  does  not  appear  likely  that  a  perfectly 
uniform  temperature  distribution  can  be  achieved.  The  second 
approach  is  capable  in  principle  of  athermalizing  an  infinitely 
long  rod  or  a  thin  disc.  But  in  any  rod  of  finite  length  there 
will  be  end  effects  which  depend  on  the  magnitude  of  the  thermal 
gradient.  A  method  for  calculating  these  end  effects  is  developed 
in  Section  6. 

Relatively  uniform  temperature  distribution  is  desirable 
for  smother  reason.  Both  temperature  distribution  and  distribution 
of  inversion  are  governed  by  the  distribution  of  absorbed  pumping 
energy.  If  the  latter  is  uniform,  both  of  the  former  will  be 
uniform  or  very  nearly  so.  And  uniformity  of  inversion  is 
important  to  give  a  laser  emission  wavefront  which  is  uniform 
in  intensity  as  well  as  in  phase .  Uniformity  in  intensity  is 
desirable  to  give  the  best  possible  far-field  pattern,  and  to 
reduce  non-linear  effects  in  the  laser  material  which  in  turn 
will  cause  phase  distortion  o.:  the  wavefront.  The  import  of 
non-linear  effects  is  discussed  in  Section  5. 

The  calculation  of  temperature  profiles  in  Section  4 
was  originally  intended  to  (a)  develop  the  method  of  calculation, 

(b)  show  the  magnitude  of  the  temperature  gradient  problem,  and 

(c)  provide  input  data  for  the  end  effect  calculations.  The 
calculations  showed  that  an  unclad  rod  of  radius  R  has  a  minimum 
temperature  at  or  near  the  surface,  with  a  substantial  gradient 
in  the  vicinity  of  R/n,  where  n  is  the  index  of  refraction 
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Water  or  glass  clad  rods  have  a  maximum  temperature  just  inside 
the  surface  (at  0.8  tt  and  0.9  R  respectively).  The  most  uniform 
temperature  distribution  is  for  the  glass  clad  rod.  An  uroxpected 
result  of  this  woak.  was  that  it  pointed  out  a  potential  method  of 
improving  substantially  the  pumping  uniformity .  That  method  is 
to  use  a  cladding  of  about  the  same  index  as  that  of  the  laser 
core  glass  bet  of  a  thickness  somewhat  less  than  that  chosen  for 
the  example  (0.5  R) .  Additional  calculations  would  be  necessary 
ho  establish  the  optimum  cladding  thickness  and  index. 

In  Section  5  the  formalism  is  developed  to  allow  calcu¬ 
lation  of  index  of  refraction  changes  due  to  electrostriction, 
photoelashic  effect,  and  Kerr  effect  when  the  laser  field 
distribution  in  the  rod  is  given. 

End-effecc  studies  were  undertaken  to  establish  under 
what  conditions  they  are  negligible,  and  what  form  they  take  when 
they  are  non-negiigible .  In  a  typical  rod  configuration,  the  end 
effects  become  important  only  when  the  glass  has  been  athermalized . 
However,  in  this  important  case  the  end  effects  can  result  in 
several  wavelengths  distortion  of  the  laser  wavefront.  For  the 
unclad  rod,  the  steep  temperature  gradient  gives  rise  to  an  optical 
path  distortion  which  is  not  proportional  to  the  temperature  distri 
bution.  Therefore,  it  does  not  lend  itself  to  correction  by 
adjusting  the  glass  parameters.  However,  in  cases  where  the 
temperature  gradient  is  less  severe,  the  calculations  indicate 
that  it  would  be  possible  to  compensate  for  the  end  effects 
almost  completely  by  proper  choice  of  glass  parameters . 
Specifically,  a  value  of  is  chosen  which  is  slightly  less 
negative  than  the  value  which  is  optimum  for  an  infinitely 
long  rod.  The  difference,  c^-a^  ,  depends,  of  course,  on  the 
length  of  the  rod. 

Analysis  of  the  limiting  cases  of  infinitely  long  rods 
(plane  strrin)  and  thin  discs  (plane  stress)  provides  insight  into 
the  athermalization  problem.  The  analyses  of  Quelle1  and  of 
Snitzer3  have  been  utilized  to  compute  the  thermal  coefficient 
of  optical  path  for  several  laser  glasses  in  rod  and  immersed 
disc  cavities.  This  work  confirms  that  athermalization  is 
possible  for  a  long  rod,  but  that  it  is  substantially  easier 
to  accomplish  for  an  immersed  disc .  The  required  indices  for 
the  immersion  fluid  are  close  to  that  of  water.  Methods  for 
removing  the  stress-birefringence  effect  are  given. 

In  all  phases  of  the  work  it  was  assumed  that  essentially 
a  single  pumping  pulse  is  to  be  used.  That  is,  no  considerations 
of  cooling  time,  or  steady  state  temperature  distribution  were 
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entertained.  enaiyses  would  therefore  apply  to  cases  in  which 

the  temperature  of  the  rod  returns  to  ambient  before  the  next  shot 
and  in  which  the  absorption  of  pumping  light  provides  the  only 
heating  effect  during  the  time  period  of  interest,  namely  the 
period  of  lasing. 
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2.  REVIEW  OF  PREVIOUS  WORK 


The  problem  of  calculating  optical  path  differences 
induced  by  temperature  gradients  in  pumped  laser  rods  has  been 
formulated  by  Snitzer  ,3  The  analysis  indicated  that  there  exist 
relations  between  the  thermal  expansion  coefficient  a,  the  index 
coefficient  c^,  ar.d  other  glare  constants  which  will  yield  an 
athermalized  condition.  In  the  writing  of  a  proposal4  for  the 
present  contract  this  formulation  was  improves  and  extended  to 
include  both  radial  and  tangential  polarization. 

Walsh6  has  studied  the  photoelastic  effects  in  a  long 
laser  rod  under  the  conditions  of  constant  heat  production  and 
surface  cooling.  E->  has  also  estimated  the  distortion  of  the 
end  face  under  these  conditions.  Quelle6  has  formulated  the 
calculation  of  optical  path  differences  and  hets  shown  mathe¬ 
matically  that  the  athermalization  condition  can  be  independent 
of  the  temperature  distribution  if  a  Pockels  type  glass  is 
assumed.  This  glass  has  a  zero  value  of  the  stress-optical 
coefficient . 
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The  distribution  of  pumping  energy  in  a  laser  rod  has 
been  considered  by  several  authors .  McKenna7  treated  the  case 
of  an  absorbing  rod  with  a  transparent  sheath,  the  absorption 
coefficient  having  a  given  value  for  all  pumping  wavelengths. 

Sooy  and  Stitch®  also  tormulated  the  calculation  of  the  pumping 
energy  distribution  in  a  laser  d  with  three-dimensional  pumping. 
Skinner®  studied  experimentally  the  distribution  of  pumping 
energy  in  a  ruby  rod  using  two-dimensional  pumping.  Borrelli 
and  Charters10  calculated  pumping  energy  distribution  in  unclad 
glass  laser  rods . 

Experimental  determination  of  thermal  distortion  of 
wavefronts  has  been  accomplished  by  Welling,  Bick  rt,and  Andresen.11 
They  have  shown  that  single  flashlamp  pumping  produces  large 
asymmetries  in  the  optical  path,  but  that  four  flashlamp  pumping 
produces  a  nearly  circular  symmetric  pattern.  They  have  3hown 
experimentally  the  development  of  the  path  differences  during 
the  pumping  pulse,  and  their  dependence  on  diameter,  doping 
concentration,  and  surface  finish  for  ruby  and  neodymium  glass. 

Work  at  American  Optical  Company  under  contract 
Nonr  3833(00)  has  been  directed  in  part  toward  the  attainment  of 
a  high  tical  quality  glas  -  which  will  also  have  a  low  thermal 
distortion  coefficient.  Data  for  a  few  of  the  glasses  studied 
under  this  contract  are  included  in  Section  3* 
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3.  THERMALLY  STABLE  CAVITIES 


A  serious  problem  in  large  lasers  is  the  thermal 
distortion  of  the  cavity  due  to  non-uniform  pumping.  Glasses 
can  be  made  with  an  index  of  refraction  variation  of  no  greater 
than  ±  0.8  x  10"6  across  a  25  mm  section,  which  corresponds  to 
2.5  fringes  ^.n  a  one  meter  length  rod.  However,  when  the  glass 
is  pumped,  substantial  changes  in  index  of  refraction  occur.11 
The  result  is  the  frequently  encountered  frustration  in  attempting 
to  produce  a  high  intensity,  diffraction  limited  beam.  Just  above 
threrhold  narrow  beams  are  possible,  but  for  increased  pumping 
the  thermal  distortions  increase .  More  energy  is  produced  but 
in  a  wider  beam.  If  a  strong  mode  selection  scheme1 3  is  used 
the  output  rapidly  levels  off  with  increased  pumping. 

In  order  obtain  diffraction  limited  beams  from  the 
end  of  laser  devices  the  cavity  should  not  distort  while  puxnilng. 
If  the  temperature  of  the  laser  rod  changes  uniformly  during 
pumping  there  will  not  be  any  distortion  of  the  wavefront  within 
the  rod.  By  placing  flashtubes  symmetrically  about  the  red  a  non- 
uniform  azimuthal  heating  can  be  avoided.  A  non-uniform  radial 
heating  of  the  rod  is  almost  inevitable  i f  reasonable  rod  sizes 
and  neodymium  concentrations  are  used.  However,  by  the  proper 
design  of  the  total  resonant  cavity  and  with  properly  chosen 
glass  compositions  one  can  minimize  the  effect  of  non-uniform 
heating . 


Due  to  pumping  of  the  laser  a  temperature  gradient  is 
produced  from  the  center  to  the  edge .  The  distortion  of  the  cavity 
arises  because  of  three  factors .  The  change  in  temperature  leads 
to  an  elongation  because  of  the  finite  coefficient  of  expansion  of 
the  glass.  With  a  change  in  temperature  the  index  of  refraction 
in  general  also  changes.  Finally,  thermal  gradients  wiunin  the 
glass  produce  stresses  which  result  in  both  a  change  in  index  of 
refraction  and  birefringence.  The  expansion  coefficient  a  and 
he  thermal  coefficient  of  index  are  functions  only  of  t\a 
amperature  of  the  point  under  consideration,  but  the  indices  of 
refraction  produced  by  the  strains  are  dependent  on  the  shape  of 
the  rod  and  the  details  of  the  temperature  distribution. 

For  glasses,  the  index  of  refraction  always  increases 
in  compression,  but  can  be  positive  or  negative.  The  effects 
producing  the  index  change  can  be  made  to  cancel  each  other  by 
an  appropriate  glass  composition  with  the  correct  negative  value 
of  c^. 


There  are  two  geometries  with  cylindrical  symmetry 
that  are  of  interest  for  which  approximate  solutions  can  readily 
calculated.  One  is  the  long  rod  in  which  the  length  L  is  much 
larger  than  the  radius  a.  The  other  is  the  disc  configuration 
with  L  «  a.  The  description  is  similar  to  that  gwen  by 
F.  W.  Quelle.1 

For  the  solid  rod,  the  total  optical  pathlength  Pr(r) 
for  a  typical  ray  parallel  to  the  axis  and  displaced  a  distance  r 
from  the  center  and  with  its  plane  of  polarization  in  the  radial 
direction  is  given  by 

\-<r>  -  nL  { 1  +  [v  -(s  <ez  +  ee>  +  7  er)  j  }  '  (1> 

where  T  is  che  difference  in  temperature  between  the  center  and  a 
point  at  a  distance  r  from  the  center.  A  cylindrical  coordinate 
system  is  used  with  the  z  direction  along  the  rod  axis.  The  e’ s 
are  the  strains.  The  quantities  q/v  and  p/v  are  the  strain-optic 
coefficients  which  relate  the  change  in  index  of  refraction  to  the 
strains  in  the  directions  parallel  and  perpendicular,  respectively 
to  the  plane  of  polarization  of  the  light.  The  notation  is  the 
same  as  used  by  Morey.13  The  corresponding  expression  for 
tangential  polarization  is 

>e(r)  -  nL  1 1  +  |anT  -  (£  +  er)  +  |  e0)j  j  .  (2) 

For  an  isotropic  medium  the  principal  strains  are 
related  to  the  stresses  a  by14 

er-E-i  [or  -  s(a9  +  <g] 

e0-E"1  [°e  *  s(°r  +  °z^]  <3> 

ez.E-x  [<rz  -  s(or  +  cre)  ] 

where  £  is  Young’s  modulus  and  s  Poisson’s  ratio. 


For  a  long  rod  with  the  end  effects  neglected,  the 
problem  is  one  of  plane  strain.  If  the  ends  are  free  of  traction 
so  that  they  can  move  in  response  to  the  heating,  tho  condition 
that  applies  is1 6 


a 

z 


a  +  oa. 
r  9 


(4) 
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In  arriving  at  Eqs .  (7)  -  (9)  for  the  long  rod  end 
effects  were  neglected.  If  the  rod  is  very  long  the  end  effects 
are  small  and  the  problem  reasonably  approximates  that  of  plane 
strain.  On  the  other  hand,  if  the  length  is  much  shorter  than 
the  diameter  the  "end  effects"  predominate  and  the  constraints 
correspond  fairly  closely  to  plane  stress. 

Analysis  of  these  equations  shows  that  there  are  three 
ways  to  achieve  simultaneously  athermalization  and  freedom  from 
stress-birefringence .  They  are: 

(1)  Athermalize  for  either  the  radial  or  tangential 
polarization.  This  would  also  require  mode  selection  so  that 
only  one  polarization  is  able  to  oscillate . 

Conditions  Long  Rod  D*sc 

Radial  APr(r)  *  0  (Eq.  7)  APr'(r) 

Polarization 

Tangential  APg(r)  ■  0  (Eq.  8)  APg'fr) 

Polarization 

(2)  Obtain  a  Pockels  glass,  i.e.,  a  glass  for  which 

p  =  q.  Quelle1  has  shown  that  the  condition  for  athermalization 
for  a  long  rod  is  then 

n{an  +  Sr(1-2s)  v}’°  (rod)  (16) 

This  esqpression  is  obtained  from  QuelleT  s  equation  by  noting  that 

B±  m  jj.  [(1-s)  7  -  s  -  £  (l-2s)  ?.  for  a  Pockels  glass. 

Equation  (16)  can  also  be  derived  from  Eq.  (9)  by  letting  p  ■  q. 

For  the  stacked  disc  case,  athermalization  of  a  Pockels 
glass  dictates  the  following  relation 

“r.  +  °  I2?”  *1+S)  +  v  }  “  0  (disc)  (17) 

(3)  Rectify  the  stress -birefringence  effect. 

On  passing  once  through  a  rod  with  a  radial  temperature 
distribution  a  phase  difference  [APr(r)  -  APg(r)]  is  introduced 
between  the  radially  polarized  ray  and  the  tangentinlly  polarized 
ray  at  the  same  ra'  .us  r.  Now  if  the  two  rays  pass  through  a  90° 


-  0  (Eq.  12) 

(15) 

-  0  (Eq.  13) 
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polarization  rotator,  the  radially  polarized  ray  will  bo  converted 
to  tangential  polarization,  and  vice-versa.  If  the  two  rays  are 
passed  again  through  the  same  rod,  or  through  an  identical  rod,  the 
phase  difference  between  the  two  rays  will  be  removed. 

For  a  single  rod  the  cavity  is  shown  in  Fig .  1 .  On  each 
round  trip,  the  light  pauses  through  the  45°  Faraday  rotator, 
thereby  undergoing  a  90°  rotation  of  the  plane  of  polarization. 


Figure  1.  Schematic  of  a  laser  with  a  rectifying  Faraday  rotator. 


Another  method  cam  be  used  when  there  are  two  identical 
laser  rods  in  either  an  oscillator  or  am  amplifier .  It  employs  a 
90°  optical  rotator  between  the  two  laser  rods,  as  shown  in  Fig.  2. 


Mirror  90°  Mirror 

Rotator 

Figure  2.  Use  of  a  90°  optically-active  rotator  to  rectify 
stress-birefringence  in  a  laser  system. 


This  rotator  can  be  of  crystalline  quartz,  for  example,  cut 
perpendicular  to  the  optic  atxis.  Other  optically  active 
materials  can  be  used. 
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In  these  cases ,  the  birefringence  is  automatically 
cancelled,  and  the  athermalization  conditions  are  found  from 
Eqs.  (9)  and  (14): 


rod 

disc 


(18) 

(19) 


The  possibilities  of  achieving  athermalization  by  these 
approaches  are  explored  more  fully  below .  At  this  point  some 
general  comments  are  in  order. 


The  conditions  for  athermalization,  Eqs .  (15)  -  (19)> 
can  be  regarded  as  requirements  on  the  ratio  c^/a.  This  is 
because  s,  p/v,  and  q/v  do  not  depend  strongly  on  the  glass 
composition.  Bridgman18  measured  several  glasses,  including 
lead  silicates,  bcrosilicates,  and  alkali  earth  silicates  and 
found  that  Poisson’s  ratios,  varies  from  0.19  to  0.26.  Furthermore, 
the  quantities  in  braces  in  Eqs.  (7)  -  (9)>  (12)  -  (14)  are  not 
sensitively  dependent  on  s.  A  nurtiber  of  references  on  the  strain- 
optic  coefficients  are  given  in  Morey.14  Pockel*  s  measurements 
indicate  p/v  »  q/v  =  0.42  for  a  silicate  glass  containing 
approximately  75$  by  weight  of  PbO.17  For  the  more  common 
glasses  the  strain-optic  coefficients  are  smaller*  q/v  decreases 
more  rapidly  than  p/v.  A  light  flint  silicate  (54.3  wt  $  SiOa, 

33  PbO,  1.5  Ba03,  3  NaaO,  8  KjO)  has  the  values  p/v  -  C.306 
and  q/v  =  0.213,  and  for  a  borosilicate  crown  (68.2  wt  $  SiOa, 

10  Ba03,  10  NaaO,  9.5  KgO,  2  AlaOa)  p/v  =  0.269  and  q/v  -  0.147.19 

To  obtain  an  estimate  of  the  required  on  to  reduce  the 
various  AP* s  to  zero,  typical  values  for  the  parameters  are  assumed. 
The  strain-optic  parameters  are  taken  as  p/v  =*0.3  and  q/v  »  0.21 
and  Poisson’s  ratio  as  s  =  0.25.  A  linear  expansion  coefficient 
of  a  =  10~6/°C  is  assumed.  If  the  temperature  T  varies  as  rp, 
then  R/T  =  l/(p  +  2);  for  a  quadratic  dependence  of  T  on  the  radius, 
p  =  2  and  R/T  =1/4.  With  these  values  for  the  parameters  the 
braces  in  Eqs.  (7)  -  (9)  are  equal  to  zero  for  0^  =  -42  x  10"7/°C, 
-34  x  1C~7 /°C,  -38  x  10“7/°C,  respectively.  For  the  stacked  discs 
in  air,  the  braces  in  Eqs.  (12)  -  (14)  are  equal  to  zero  for 
on  =  -56  x  10-7/°C,  -50  x  10~7/°C,  -53  x  10-7/°C,  respectively. 

Note  that  the  required  values  of  0^  are  all  negative.  In  most 
silicate  glasses  is  not  sufficiently  negative  to  satisfy  any 
of  the  conditions  for  zero  pathlength  difference  between  rays  at 
the  center  and  at  the  edge  of  the  rod .  If  the  temperature  is 
higher  at  the  edge  than  at  the  center  the  pathlength  is  greater 
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at  the  edge  and  furthermore,  radially  polarized  light  propagates 
more  slowly  than  light  with  tangential  polarization.  If  the 
quantity  in  the  braces  is  not  zero  but  is  reduced  to  10"7/oC, 
three  fringes  at  1 .06  m-  would  appear  between  center  and  edge  if 
T  ■  10° C  for  a  meter  length  rod  made  from  a  material  of  index 
n  -  1.5. 


At  low  temperature  On  is  negative  because  expansion  of 
the  glass  reduces  the  amount  of  polarizable  material  per  unit  volume-, 
at  high  temperatures,  this  is  more  than  offset  by  the  increased 
thermal  population  of  higher  vibration  states  of  the  ground 
electronic  state,  which  appears  to  shift  the  fundamental  ab¬ 
sorption  band  of  the  material  to  longer  wavelengths  and  thereby 
increases  the  index  of  refraction.  The  temperature  T0  at  which 
an  ■  0  is  a  function  of  the  glass  composition  and  wavelength. 

It  increases  with  increasing  k  through  the  visible  and  near 
infrared.  For  most  glasses  T0  is  less  than  room  temperature 
at  1  |x,  but  it  is  above  500°  K  for  a  few,  such  as  the  silicates 
with  high  barium  content. 

Molby20  and  Prod’ homme21  give  values  of  an  in  the  visible 
region  of  the  spectrum  for  various  glasses.  Measurements  were  made 
in  this  laboratory  at  1 .06  p.  by  using  a  thin  plate  of  the  glass  to 
be  measured  as  the  end  reflector  of  the  laser  cavity. 22  Laser 
emission  occurs  at  those  wavelengths  at  which  the  plate  thickness 
and  index  give  an  odd  number  of  k/4.  From  the  shift  in  wavelength 
of  the  laser  lines  as  a  function  of  temperature  of  the  plate  and 
by  independent  measurements  of  the  expansion  coefficient  of  the 
glass,  cn  was  determined.  Room  temperature  values  for  various 
glasses  were  between  29  Y  10_7/°C  and  -4l  x  10_7/°C.  The  positive 
values  were  for  dense  lead  silicates.  High  barium  content  silicates 
gave  the  large  negative  values. 

DISCUSSION 

(l)  Athermalization  for  radial  or  tangential  polari¬ 
zation.  It  is  necessary  to  mode-select  a  field  distribution  with 
axially  symmetric  polarization,  such  as  the  TE0l  or  TM^  modes. 

For  example,'  the  TM01  mode  is  radially  polarized  in  the  electric 
vector . 


Its  intensity  distribution  is  [  (ujr/a)]3 ,  where  Uj  is 

the  first  zero  at  3.832  of  the  first  order  Bessel  function  .  If 
one  end  of  the  glass  rod  is  cut  in  the  form  of  a  Brewster  angle 
cone  and  a  plane  mirror  placed  at  the  proper  distance,  all  the 
modes  except  the  family  TMom(m  si)  are  discriminated  against. 

The  lowest  order  member  of  this  family  could  then  be  selected 
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T 

by  other  means,  such  as  a  limited  aperture  at  the  other  reflector. 

Thi3  approach  has  several  disadvantages.  First,  it  is  difficult 
to  make  cones  with  good  optical  quality.  Second,  the  beam  is  a 
hollow  cone,  which  is  not  as  desirable  as  the  lowest  order  HE^ 
mode.  Third,  the  requirements  on  an  depend  on  R  and  hence  on  the 
details  of  the  temperature  distribution. 

(2)  For  a  Pockels  gl«.ss,  the  required  values  of  an  are 

calculated  as  follows.  In  Eqs.  (16)  and  (17)  substitute  the 
typical  values  p/v  =  0.45  and  s  «  .25  >  n  ■  1.5>  “  1*33* 

From  Eq.  (16):  an  -  -0.6  a  (rod) 

From  Eq.  (17):  =  -.367  a  (disc) 

The  original  Pockels  glass  has  a  positive  value  of  c^.  Efforts 
to  obtain  a  negative  value  for  an  in  a  Pockels  glass  will  be 
reported  under  contract  Nonr  3835(00). 

(3)  Rectification  of  the  stress-birefringence  effect 
appears  at  present  to  be  the  most  promising  approach.  In  order  to 
compare  the  various  possible  configurations  and  available  glasses, 
the  thermal  coefficient  of  optical  pathlength,  _1_  dP  ,  is 

nL  dT 

tabulated  below.  The  expressions  are  obtained  from  Eqs.  (9)  and 
(14)  for  the  average  optical  path  P,  for  the  rod  and  immersed  disc 
cases  respectively. 

r  -  I 

i 

1 
| 
| 


M  an  +  a  .125  +  .1162 J 


1  dP 


(rod)  -JrS 


a  + 


nL  dT  n  2(l-s) 


[(3-5s)  2.  +  (1-3.)  J] 
»  a  +  .383  a 

n 

/  .  x  1  '  o,  A  „  ln~n'  fi  R  +  ill  ll 

<disc>  ^ii"an+aL~(1’8)  +“v  +  ~vj 


JL 


I 


T 
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TABLE  I 


Thermal  Coefficient  of  Optical  Path  for  a  Long  Rod,  a  Disc  in  Air 
(n'  ■  1.0),  and  a  Disc  Immersed  in  Water  (n#  ■  1.55) 


(Rod) 

(Disc) 

Glass 

n 

a 

a 

n 

_L  dP 

nL  dT 

_1_  dP' 
nL  dT 

n'-l.O  n  '**1 .33 

5855# 

1.51 

105‘  10- 

7  -22*  10“7 

11.8*10- 

7  14*  10“  7  -3 . 1  •  10"  7 

1205 

1.52 

122 

-28 

18.8 

58.2  5.5 

1204 

1.57 

119 

-41 

4.6 

26.8  -4.4 

1265 

1.52 

112 

-52 

11 

28.8  -1.5 

1270 

1.55 

114 

-56 

7.7 

26.3  -4.2 

1276 

1.62 

111 

-28 

14.6 

Vrf 

00 

VO 

Co 

Constants  used:  £  B 

V 

.30,  i  -  .21, 

s  =  .25, 

except  3835. 

♦For  3835  £ 

V 

-  .25, 

-  -  .15,  8  = 

V 

.225. 

Several  conclusions  are  apparent  from  this  table.  For 
the  rod  case  a  number  of  experimental  glasses  have  a  lower  thermal 
coefficient  than  the  standard  5835  glass,  number  1204  being  a 
factor  of  three  lower.  lfce  disc  in  air  is  considerably  poorer 
than  the  long  rod.  But  the  disc  immersed  in  water  (n'  «  1.55) 
appears  very  attractive.  The  value  -1.5  *  10"7  for  glass 
number  1265  means  that  for  aim  length  of  glass,  nL  =  1.52  m, 
a  20° C  temperature  differential  between  center  and  edge  of  the 
discs  would  produce  an  optical  path  difference  of 


AP'  -  nL  dP'AT  -  1.52  •  103  (-1.5  '  10-7) 
dT 

•  ^  *  4.56 


10" 4  cm  «=  4.56  nm 


Furthermore,  the  numbers  in  the  last  column  of  the  table 
suggest  that  it  may  be  possible  to  find  either  (a)  a  glass  compo¬ 
sition  which  is  precisely  athermalized  /  .  — ,  \ 

(  _L  dP  =  0  ) 

\  nL  dT  / 

for  the  water  immersed  disc,  or  (b)  an  immersion  liquid  with 
index  n'  which  will  exactly  athermalize  an  existing  glass. 
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For  the  constants  p/v,  q/v,  and  s  assumed  above,  the 
required  index  nL  for  the  immersion  liquid  is  given  by 

+  .116  a  I 

1.25  a  J  ' 

The  required  indices  for  a  few  glasses  are  listed  below 
Glass  Required  n ' 

5835  1.38 

1203  1 .38 

1263  1 .34 

The  required  indices  are  all  within  the  range  attainable 
with  water  solutions  of  inorganic  materials. 

The  general  expression  of  the  required  index,  n ' ,  is 
obtained  from  Eq.  ( 19) 

(Disc)  n'  -  [-1  +  1+s  +  E  +  1^5.  S.1  (20) 

1«  [a  2  v  2  vj 

The  assumption  has  been  made  that  the  index  of  the 
liquid  does  not  change  during  the  pumping  cycle.  Since  there 
can  be  very  little  heat  transfer  to  the  liquid  from  the  laser 
glass  during  the  pumping  period,  the  assumption  implies  that  the 
liquid  is  non-absorbing  to  the  pumping  radiation . 
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4.  CALCULATION  OF  TEMPERATURE  PROFILE 


4.1  ASSUMPTIONS  AND  DEFINITIONS 

If  we  may  assume  that  the  temperature  change  in  a  laser 
rod  is  proportional  to  the  radiation -heating  effect  due  to  ab¬ 
sorption  of  radiation  from  an  adjacent  flashtube,  then  the 
calculation  of  the  temperature  profile  reduces  to  a  calculation 
of  the  distribution  of  heat  production  in  the  laser  rod. 
Measurements  on  a  large  variety  of  glasses33* 84  show  that  the 
specific  heat  capacity  is  sufficiently  constant  over  the  tempera¬ 
ture  range  expected  in  the  laser  rod  so  that  the  proportionality 
may  be  assumed. 

A  computer  program  has  been  written  by  which  to  calculate 
the  distribution  of  heat  production.  The  program,  written  in  the 
FORTRAN  language  for  an  IRM-36O-3O  computer,  allows  for  three 
different  homogeneous  media  to  be  placed  between  the  flash tube 
and  the  laser  rod.  Cylindrical  geometry  and  polished  refracting 
surfaces  are  assumed  so  that  all  boundaries  between  media  form 
infinitely  ’ong,  cylindrical  surfaces  concentric  with  an 
infinitely  long,  cylindrical  laser  rod. 

It  is  assumed  that  the  flashtube  and  reflector  produce 
a  volume  of  isotropic  radiation  (called  the  source  hereafter)  in 
which  the  rod  and  its  cladding  media  are  immersed  (no  variation 
with  azimuth  about  the  axis  of  the  laser  rod) . 

It  is  assumed  that  the  radiant  flux  absorbed  by  the 
laser  material  goes  entirely  into  exciting  the  doping  ions  to 
higher  energy  levels  plus  a  radiative  decay  that  produces  the 
laser  radiation. 

It  is  assumed  that  each  photon  absorbed  results  in  a 
photon  of  laser  radiation  while  the  rest  of  the  energy  is  converted 
to  heat. 


It  is  assumed  that  the  transmissivity  of  the  laser 
material  does  nor  change  significantly  with  change  in  population 
density  of  excited  states  of  any  of  the  ions . 

The  effect  of  multiple  passes  of  the  radiation  through 
the  rod  are  not  taken  into  account  by  the  program.  This  effect 
will  be  most  noticeable  for  weakly  absorbing  and/or  small  diameter 
rods.  Partial  compensation  for  this  emission  can  be  carried  out 
by  (1)  increasing  the  overall  energy  content  of  the  flashtube 
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spectrum  and  (2)  adjusting  the  spectrum  to  favor  the  weakly 
absorbed  wavelengths. 


Parameters  which  affect  the  distribution  of  heat 
production,  and  which  must  be  specified  at  the  time  the  program 
is  run,  are:  (1)  the  product  of  the  radius  times  the  absorption 
coefficient  for  the  laser  rod  and  each  layer  of  its  cladding, 

(2)  the  refractive  indices  of  the  media,  and  (?)  the  spectral 
radiance  of  the  source . 


The  following  definitions  have  been  adapted  from 
references  (25,  26): 

radiant  flux  -  rate  of  transfer  of  radiant  enercv,  specified  in 
watts 

radiant  transmittance  -  or  simply  transmittance,  is  the  ratio  Wt/»i 
of  the  transmitted  flux  tc  the  incident  radiant  flux. 


spectral  transmittance  -  the  radiant  transmittance  evaluated  for 
a  particular  wavelength  of  the  incident  energy. 

transmissivity  -  the  value  of  the  internal  transmittance  for  a 
unit  thickness  of  a  non-diffusing  substance. 

absorption  coefficient  -  Bouger'  s  law  states  that  equal  layers  of 

an  imperfectly  transparent  material  will  absorb  equal  fractions 
of  the  radiant  energy  entering  them.  If  this  law  is  applied 
to  infinitesimal  layers  and  if  the  losses  are  integrated,  the 
radiant  flux  at  a  distance  x  within  the  medium,  for  the 
wavelength  is  found  to  be  P  =  PQ  exp(-a(k)x),  where  PQ 
is  the  flux  of  wavelength  k  entering  the  first  layer  and 
a(M  is  the  absorption  coefficient. 

radiance  -  the  radiant  flux  or  power  per  unit  solid-angle-in-the- 
direction-of-a-ray  per  unit  pro jected-area-perpendicular-to- 
the-ray.  It  has  the  same  value  at  any  point  along  the  ray 
within  an  isotropic  medium  in  the  absence  c£  looses  by 
absorption,  scattering,  or  reflection,  (it  may  be  helpful 
to  note  that  radiance  is  analogous  to  the  photometric 
quantity  luminance  or  photometric  brightness) . 

spectral  radiance  -  the  radiance  at  any  specific  wavelength. 

Typical  units  are  watts/cm3 -ste radian  per  millimicron. 

spectral  irradiance  -  radiant  flux  incident  per  unit  area  of  a 
surface  at  any  specific  wavelength.  Typical  units  are 
watts/cm3  per  millimicron. 
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4.2 


DERIVATION  OF  THE  BASIC  EQUATION 


tical  pumping  of  laser  rods  subjects  the  laser  material 
to  a  distx j.-.ution  of  energy  density;  power  is  extracted  from  the 
radiation  field  in  the  form  of  heat  and  laser  radiation.  Consider, 
then,  the  simple  one -dimensional  case  of  an  elemental  volume  of 
cross  section  dA  and  length  d i.  Let  the  energy  density  due  to 
radiation  in  the  wavelength  interval  dX  centered  on  the  wavelength 
>.  be  uniform  and  of  strength  dE  in  the  volume  dA  •  di.  Let  it 
arxse  from  collimated  radiation  propagating  along  the  direction 
di.  Then  the  radiation  in  the  volume  dA  •  di  will  move  a 
distance  At  •  c/nl  in  time  At,  where  nx  is  the  refractive  index 
of  the  medium  a.id  c  is  the  velocity  of  light  in  vacuum.  If  the 
medium  is  absorbing  with  a  linear  absorpu.  :i:  coefficient  a x(X), 
the  energv  W  -  dE  •  dA  •  di  will  be  diminished  to 

W'  -  W  j,<xx  cAt/nx  (21) 

in  its  passage  through  the  medium.  The  energy  absorbed  will  be 

w  .1  W-W'  =  W  ,  (22) 

or,  for  a  short  time  intervi.  . 

w  a  wt,  —  At.  (23) 

-  Uj 

That  is,  the  energy  absorbed  per  unit  of  time,  or  the  r  *r  loss 
will  be 

v/At  ■  Wa1c/n1 ,  (24) 

Ibis  power  loss  occurs  in  the  volume  dA  •  di,  so  the  power  loss 
per  unit  volume  is 

dp'  *=  w/At/(dA  •  di)  «*  c^cdE/nj  .  (23) 

Of  chis  loss,  only  the  fraction  (1-X/Xe)  goes  into  thermal 
agitation  of  the  medium,  where  Xe  is  the  wavelength  of  the  laser 
radiation.  The  heat  production  per  unit  volume  due  to  the  radi¬ 
ation  in  the  wavelength  interval  dX  is 

dp  -  ax  — -  dE(l-X/Xe)  watts/cm3  (26) 

where  dE  is  termed  the  spectral  energy  density. 
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The  spectral  energy  density  at  a  point  M  can  be  found 
as  follows:  Let  the  point  be  contained  in  a  surface  element  dA 
through  which  the  spectral  irradiance  is  HM(7v)  .  The  partial 
energy  density  due  to  that  radiation  is 

dE(M  =  lyx)  n!  (x)/c  (27) 

The  total  spectral  energy  density  for  incoherent  radiation  would 
be  found  by  adding  up  all  such  contributions  at  the  point  M. 
Suppose,  however,  that  the  spectral  irradiance  Kfo(X)  is  due  to 
radiation  contained  in  a  solid  angle  do)  centered  around  the  normal 
a  to  the  surface  element  dA.  Then  the  spectral  irradiance  Hm(^) 
at  M,  being  the  flux  dP  passing  through  dA  could  be  expressed  as 

H^k)  =  dF/dA  -  NM(>.,o)  do,  (28) 

where  NM(k,c)  is  the  spectral  irradiance  at  M  due  to  radiation  in 
the  wavelength  interval  dk .  That  is,  NM(k,a)  is  the  radiant  flux 
per  unit  solid  angle  in  the  direction  of  the  normal  per  unit 
projected  area  perpendicular  to  the  normal.  The  summation  of 
all  such  contributions  at  the  wavelength  k  becomes 

dE„(k)  =  f  NM(*,o)  do.  (29) 

Now,  NM(k,tr)  is  the  spectral  radiance  at  the  point  M 
within  the  laser  rod  and  is  rot  known  directly.  But,  in  the 
homogeneous  medium  of  the  rod  with  index  nx(k)  and  absorption 
coefficient  ax (k)  the  spectral  radiance  along  a  given  ray 
diminishes  according  to  the  exponential  decay  law  expf-Oji). 

Thus,  the  spectral  radiance  Nm  in  a  given  direction  at  M  must 
be  a  diminution  of  the  spectral  radiance  Nx  in  the  same  direction 
just  inside  the  surface  of  the  laser  rod; 

hm(*,o)  =  ^(*,0)  exp  (a)]  ,  (30) 

where  dx(o)  is  the  pathlength  of  the  ray  from  the  surface  to  the 
point  M. 


In  the  absence  of  reflection  los,  -s,  the  spectral  radiance 
Nj  would  be23 

Hx  -  NB(nx/na)3,  (31) 

where  Ns  is  the  spectral  radiance  in  the  direction  specified  by 
a  and  Snell's  law,  and  ns  is  the  refractive  index  of  the  material 
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surrounding  the  laser  rod.  Nj  is  diminished  from  this  value  by 
the  Fresnel  coefficient  of  transmission  T{a). 

N  =  T (a)  Na(*,fi)  e"ai(X)di('T'.  (32) 

w  n  « 

“2 

The  ray  may  be  traced  back  through  the  various  cladding  media ,  in 
this  way,  to  the  isotropic  source  volume.  Such  a  source,  by 
definition,  has  a  spectral  radiance  N0(7\)  which  is  independent 
of  the  angle  of  exit  from  the  source . 


For  a  general  jth  ray, 

N  (>,<"  )  «  T  (*)No(?0  nxz  exp  I"-  £  a  (*)d(m,j)  |  ,  (33) 

M  D  l  m  ™  J 

where  Tj(^)  is  the  product  of  the  transmission  coefficients  at 
the  interfaces  encountered  by  the  jth  ray,  am(X)  is  the  absorption 
coefficient  of  the  m  h  medium  along  the  ray,  and  d(m, j)  is  the 
pathlength  in  that  medium  for  that  ray. 


Equations  (29)  and  (33)  yield 


ni3  C 

dEM(>)  =  /  T^(^)  N0(x)  exp 


-  £am(Md(m>:j)l  dco.  (34) 

ill  —  f 


The  integration  indicated  in  Eq.  (34)  was  done  numeri¬ 
cally  by  dividing  the  angular  space  about  M  into  angular  intervals 
and  summing  the  contributions  from  all  the  intervals .  Because  of 
symmetry  we  need  integrate  over  only  one  quarter  of  vhe  space  so 
it  becomes  convenient  to  express  da>  in  terras  of  the  polar  angles  y 
and  a ,  where  y  is  measured  from  the  axis  of  the  laser  rod  and  a  is 
measured  in  a  plane  normal  to  the  axis.  Uius, 


do  =  sinydady. 


(35) 


Let  a  be  incremented  in  the  fashion 

1/2  Act  (Act)  180° -1/2  Aa,  (36) 

and  let  y  be  decremented  in  the  fashion 

90° -t/2Ay  (Ay)  1/2  Ay.  (37) 

Then  no  ray  with  polar  angles  (a,y)  lies  in  any  of  the  planes  xy, 
yz  or  xz,  (a  measured  from  the  x-axis,  y  measured  from  the  z-axis). 
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yet  the  intervals  are  regularly  spaced  about  M.  Equation  (34) 
becomes 

n  3  (> ) 

dE„  -  -  — —  £  T.(*)N0(*)sinY  AaAy 

h  c  j  J  J 

•  e3£p  Eam(>')d(ni^)  J  (38) 

The  total  boat  production  per  unit  volume  in  the  laser 
rod  can  be  found  by  substituting  Eq.  (58)  into  Eq.  (26)  and  inte¬ 
grating  over  the  wavelength  range  of  the  radiation  from  the  source 
contributing  to  the  energy  density  . 

EJf  =  y*aa  (*)(l-yxe)  ni2(M  AaAy 

*  £  T.(>v)N0(7v)  exp  I"-  £a  (k)d(m,;j)l  siny-dk.  (39) 

3  -*  L  m  J  J 

Equation  (39)  is  the  basic  equation  for  the  computation 
of  the  distribution  of  heat  production  in  the  laser  rod.  The 
integration  over  wavelength  was  done  numerically  using  Simpson’s 
rule.  The  operations  required  to  apply  the  equation  form  a  four¬ 
fold  nest  of  loops; 

1 .  Trace  the  selected  ray  from  the  point  m(x,0,0) 
in  the  laser  rod  back  to  the  source  for  a  chosen 
wavelength.  Compute  the  transmission  coefficients 
at  each  of  the  interfaces,  and  evaluate  the  contri¬ 
bution  of  the  flux  arriving  at  P  from  the  selected 
direction  (innermost  loop). 

2.  Change  the  angles  specifying  the  ray  anc.  do 
(1)  again.  Keep  repeating  these  operations  until 
all  the  flux  in  the  interval  d^  arriving  at  M  has 
been  accounted  for. 

3.  Add  the  energy  from  this  wavelength  interval 
to  the  energy  contributed  by  other  intervals 
according  to  Simpson’ s  rule  for  numerical  inte¬ 
gration.  Repeat  the  entire  loop  for  each  of  the 
wavelength  intervals  required  to  cover  the 
operative  spectrum  of  the  source . 

4.  Repeat  the  entire  calculation  for  various 
points  M(x,0,Q)  along  the  radius  of  the  laser 
rod  selected  to  best  show  the  distribution  of 
heat  production . 
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4.3  POLISHED  LASER  ROD  WITH  CLADDING  AND  ISOTROPIC  RADIATION 

WITHOUT  MULTIPLE  REFLECTIONS 

Basic  Ray-Tracing  Equations 

Let  the  laser  rod  of  index  n^  be  immersed  'n  a  medium 
of  index  ns .  Let  the  infinitely  long  rod  be  ?  cylinder  of  radius 
Rr .  Let  the  flashtube  be  approximated  by  a  uniformly  radiant  gas 
volume  separated  from  the  laser  rod  cladding  by  a  transparent 
envelope  forming  an  infinitely  long  cylindrical  cavity  concentric 
with  the  laser  rod.  The  cavity  has  radius  Ra  and  index  na .  This 
geometry  is  illustrated  in  Fig.  3>  which  includes  a  rectangular 
coordinate  system  with  z-axis  coincident  with  the  axis  of  the 
laser  rod.  We  wish  to  consider  a  typical  point  P  in  the  laser 
rod  and  trace  rays  from  it  to  the  surface  of  the  surrounding 
gas  volume . 


Figure  3*  Geometry  of  the  laser  rod,  cladding  layers  and  source. 

The  point  P  of  observation  can  be  taken  on  the  positive 
x-axis  without  loss  of  generality;  P  =  P(x,0,0).  The  xz-plane  and 
the  xy-plane  are  planes  of  symmetry.  The  total  radiation  arriving 
at  P  can  be  obtained  by  ray  tracing  in  the  sector  y  s  0,  z  s  0  only, 
as  long  as  we  neglect  reflected  rays.  An  arbitrary  point  on  the 
surface  of  the  laser  rod,  in  this  sector,  will  be  called 
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pi (xi >  Yi >  zi ) >  where  the  equation  of  the  surface  of  the  rod  is 

Fi  (xi  >  Yx  »  zi  )  "  0  “  xi2  +  Y;a  “  Ri2  •  (40) 


The  geometrical  distance  from  Px  to  P  is 

d,  o  j\x-xj2  +  y,2  +  *xa  J  * 

The  optical  direction  cosines  of  the  ray  proceeding  from 

pi  £2.  p  aro 


Pi  =  nx 


-  “ni— ;  ri 


(41) 


(42) 


The  ray  which  is  bent  alcng  this  direction  by  refraction  at  the 
surface  has  optical  direction  cosines  (p2  ,  qa ,r2  )  related  to 
(Pi  ,qi  )  by  Snell*  s  law; 


P2  “  pi  “  ?  0Xi  '»  q2  “  qi  =  5  ay!  ’  rs  "  ri  ”  ^  Bz1  ’ 

where  £  is  an  (as  yet)  undetermined  multiplier.  Such  a  ray 
proceeds  from  the  source  point  P2 (x2 ,  y2 ,  z2 )  with 


(4?) 


X1  -  X2 


Yx  “  Ys 


Z1  “  Z2 


Pa  =  n2 


"5  ^2  ~  n2  3  ')  r2  n2  (J 


(44) 


where 


* 

d2  -  £(xi  -  x2)2  +  (y!  -  y2)?  +  (zx  -  z2)aj 
Finally,  we  add  the  conditions  that 


(45) 


Pi2  +  <*x3  +  ti2  -  n3 


p22  +  qa2  +  r22  -  n2; 


(46) 


in  order  to  complete  the  system  of  equations.  We  will  solve  the 
system  in  the  sense  that  —  given  the  paroneters  x,  ^ ,  n2 ,  Rx ,  Ra  , 
find  p8 ,  q2  ,  ra  ,  da  and  cL,  as  functions  of  px ,  qx  ,  rx  . 

We  find  x1 ,  z1  by  solving  Eqs .  (40)  and  (42) 
simultaneously.  Thus, 


Xi  -  x  +  ylPl/qi 


(47) 
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. . 


so  that 


Yx 


Rx2  =  x2  +  y12(p1/q1)2  +  2  xyx(px/qx)  +  yx2 

Vi 2  [1  +  (Pi/qx)2]  +  2  X  px/qx  yx  +  Xs  -  R,2  -  0 

.  p*/*  _  _  fj  pi  a.  _  -  J2  +  v  - x2 

1  +  (pxAx)2  Li1  +  (px/'Ji)2  )  1  +  (pi/^x)2 


(48) 

(49) 

(50) 


The  negative  root  is  chosen  in  order  to  work  only  with 
the  radiation  coming  toward  p  from  the  negative  half  space  y  ^  0. 
It  is  desirable  to  use  Eq.  (50)  whenever  Jp1  j  <  |qx  |,  obtaining 
xx  from  the  relation 

xx  -  x  -  (y  -  yx)  Pj/q.  (51) 


and  obtaining  z1  from  the  relation 
*x  -  +  yiri/<h  • 


(52) 


When  jpx i  >  |qx|,  Eqs.  (40)  and  (42)  may  be  solved  to 

yield 


(q,  /Pi  )2  v  -l  r 

1  (qx/Px 52x  ) 

3  x  V  -  *3(<j,.':-1;31 

1  1  +  (qi/Pi  )s  L 

if  sgn(pn)  » 

(qx/Pi  )2  v 

U  +  K/PxH 

jgn(x  -  xx ) 

(  tex/Px)2* 

i  +  (qi/Pw  J 

|S  ^  Rx3  -  x2(q1/p1)3  1 

i  +  (qx/Pi)2 

(1  +  (qi/Px  )2  , 

+  i  +  (qi/p,  )2  J 

if  the  sign  condition  is  not  met  by  Eq .  (55). 

In  both  of  these  cases  one  obtains  yx  and  zx  from  the  equations 


Yi  =  y  -  (x  -  Xx  )  qx/px 

sx  -  (xx  -  x)  rx/px 

Next  we  find  £  in  terms  of  xx  ,  yx ,  zx  , 


££i 

3xx 


2xx ; 


111 

^Yx 


9Fx 

2y-  3^ 


o 


(55a) 

(55b) 

(56) 
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p£  =  (-  p1  ;  qs  =  2^y1  +  q3  ;  ra  «  rx 


(57) 


By  squaring  each  of  the  relations  in  Eq .  (57)  and  adding,  we 
obtain  by  Eq .  (46), 

ns3  =  n33  +  4£3(x13  +  y12)  +  4£(p1x,  +  qxy:)  (58) 

or, 

?  _  _  Pixi  +  ViVi  ^  [  (  Pi«i  +  <hYi 

2R1S  [(  2R, 3 

In  Eq.  (59),  that  root  roust  be  chosen  which  makes  £  approach  the 
physically  correct  limit  £  -  0  as  nx  -  n2 ;  i>e,,  plus  if 
(p1x1  +  q1yl)  >  0  and  negative  otherwise. 

Equations  (41),  (45),  (57)  and  (59)  allow  us  to  calcu¬ 
late  pa ,  q2  ,  r2 ,  dx  and  dg  as  functions  of  p3 ,  q1 ,  rx  ,  which  is 
the  solution  we  sought. 


Forbidden  Rays 


Equations  (50),  (55),  (54)  and  (59)  all  contain  radicals 
which  must  be  real  for  a  physically  acceptable  solution  to  the  ray¬ 
tracing.  The  radicals  of  Eqs .  (50),  (55)  and  (54)  become  imaginary 
only  if  x  exceeds  Ra ,  a  condition  which  will  not  occur .  on  the 
other  hand,  the  limitation  that  the  radical  in  Eq.  (59)  be  real 
delineates  the  angular  region  in  which  the  rays  must  lie  in  crder 
to  be  refracted  into  the  laser  rod  toward  P  rather  than  be  totally 
reflected  by  the  interface.  The  condition  is  that 

(Pixi  +  *  Ri3  (nia  -  na3).  (6°) 


An  instructive  example  is  the  case  already  investigated 
by  W.  R.  Sooy  and  M.  L.  Stitch,  wherein  n2  =  1 .  [  "Energy  Density 
Distribution  in  a  Polished  Cylinder  of  Laser  Material,"  J.  Appl. 
Phys .  ^4.,  1719  (1965).]  In  particular,  for  rays  in  the  xy-plane, 
we  may  construct  the  typical  ray  diagram  of  Fig.  4.  By  the  law 
of  sines, 


R.  sinG 

sina,  =  — - 

1  x 


(61) 


or. 


px  =  n1cosa1 

Pi  -  ±  ni 


R33  sin30  ^ 
x3 


(62) 
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Figure  4.  Typical  ray  diagram  for  rays  in  the  xy-plane  of  a 
polished  cylinder  of  laser  material . 

A  ray  in  the  xy-plane  will  be  "forbidden"  if  [ (Rj/xJsinG  ]2  >  1. 
For  small  values  of  x  the  inequality  always  can  be  satisfied  by  an 
appropriate  value  of  0.  However,  the  maximum  value  for  sin©  it 
l/nl  (critical  refraction)  so  that  the  largest  value  x  can  have 
and  still  allow  the  relation  to  -  satisfied  for  all  p1  is 
x  »  Rj/ni.  When  x  is  larger  thau  Ax/ n x  the  region  for  which 


is  a  region  vhich  contains  no  rays  headed  towards  the  point 
p(x,0,0) .  Relation  (6C)  is  a  more  general  form  of  relation  (6j) 
since  it  applies  to  the  three-dimensional  case  with  arbitrary 

nj  <  Oj  . 

Figure  5  attempts  to  display  typical  surface  areas  of 
the  rod  which  can  irradiate  the  point  P.  Figure  5a  is  sketched 
for  x  <  Rjng/nj ,  while  Fig.  5b  is  sketched  for  x  >  . 
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point  P. 

Transmission  Coefficients 

We  follow  the  notation  of  Born  and  Wolf,  "Principles 
of  Optics,"  and  use  the  average  coefficients  for  unpolarized 
radiation.  Ihus,  from  Bom  and  Wolf,  ®  1.5*3i  Eqs.  (6l),  (60) 
and  (51)  in  that  order. 


T  -  1-R;  R  *  1/2  (R„  +  Rx) 


(64) 


tana(ei-  et)  J  R  _  ain3^-  8t) 
tana(0j+  @t)  sina(0£+  0*.) 


(65) 
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(66) 


So, 


*CO8a(0j+  0^) 

’sin(0i“  0t)"ja 

COSa(0i-  ©t) 

sin(0i+  0t) 

rfco80icos0t-  ain0isin0t}a+  {cos0icos9t+  sin0isin0t}8"| 
I  [co80£co80^+  sin0^ain0^)a  j 


‘ 

•  - 
1 


‘sin0icos0t-  cos0isin0tla 


8in0^cos0t+  cos0£sin0tj 
But  by  Snell's  law,  n^sin0^  »  ntsin0t,  so 


(67^ 


R  * 


COS 


*0t  [l-(nt/ii)asina0t]  +  sina0t  •  (nt/ni)asina0t 

i - : - : — r— - - - 


[cosQ^  \1-  (n^/ni)asina0t  +  sin0t  •  (nt/ni)sin0t]a 


cos0t  •  (nt/ni)sin0t  -  sin0t^l-  (nt/ii)a8ina0t 

cos0t  •  (nt/^i)sin0t  +  sin0tV"l  (nt/hi)asina0t 


(68) 


( 

cosa0t  [1-  (nt/ni)asina0t]  +  (nt/n£)asin40t 

T.l  | 

£cos0t  yj  1*“  (nt/hi)asina0t  +  (nt/hi)sina0t] 

■  3 

-^1-  (nt/hi)asina0t  -  (nt/n^)cos0t 

|  (69! 

V 1-  (nt/hi)asina0t  +  (nt/hi)cos0t 

A  The  angle  @t  in  Eq.  (69)  is  the  ar.gle  between  the 
normal  N  to  the  (cylindrical)  surface  between  medium  ni  and 
medium  nt  and  the  vector  al^ng  the  ray  with  direction  cosines 
(ptAt>  xt/nt).  That  is. 


A  A 

_  -  x  N  *  Pt 

CO80t  «  COS(N,  Pt)  -  — -  . 

l»llPt! 


(70) 
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For  example,  by  Eq.  (40)  for  the  surface  of  the  laser  rod,  the 
angle  o*  refraction  for  the  ray  passing  from  (xj,  ylf  zt)  to 
P(x,0,0)  is  such  that 


(71) 


2x1p1/n1  +  2y1g1/n1  ^  xxpx  +  yxgx 
2^xxa+  y1a  R1n1 


(72) 


The  angle  <f>1  must  lie  in  the  range  0  n  4>x  *  t/2,  so  we  may  take 

<t>1  «  arccos  ||(x1p1  +  y^/in^i )  \  }  (73) 

without  ambiguity  due  to  positive  or  negative  arguments. 

Similarly,  the  angle  of  refraction  #8  into  the  uodium  n8 
from  a  surrounding  medium  must  be 


08  «  arccos  {|(xaPa  +  Ya<5a)/(naRa)l  }  >  (74) 

^here  xa,  y8,  p8  and  q8  are  determined  by  the  equations  in 
section  4.2,  Basic  Ray-Tracing  Equations. 
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Flow  Diagram  of  the  Co.  .niter  Program  (polished  rod) 


IPptWSW  »tHWItlllWWt«lll*l«w*i® 


4.4 


spectral  radiance  of  the  source 


"Isotropic  radiation'1  is  realized  when  the  spectral 
radiance  N0  is  the  same  for  all  points  on  the  surface  and  for 
all  directions  of  emission  from  the  surface  of  the  source  (gas 
volume  plus  transparent  envelope  of  Fig.  3)-  Su*;h  a  source  is 
known  as  a  Lambertian  radiator  because  the  flux  emitted  per  unit 
solid  angle  (radiant  intensity)  must  be  proportional  to  the  cosine 
of  the  angle  measured  from  the  normal  to  the  surface  in  order  that 
the  raiiance  be  independent  of  the  viewing  angle. 


Goncz  and  Newell37  have  made  measurements  or  xenon 
discharge  tubes  using  an  EG&G  model  585  spectrorad’ometer.  In 
particular,  for  an  FX-47A  flashtube  (16.5  cm  arc  length,  1.3  cm 
bore,  0.4  atm  Xe  fill  pressure)  they  report  the  ratio  G  of  the 
total  light  energy  emitted  per  millimicron  by  the  lamp  to  the 
total  electrical  energy  supplied  to  the  lamp  per  discharge .  Two 
different  curves  of  spectral  output  are  given  corresponding  to 
two  current  densities  (1700  and  5300  amp/cma)  in  the  flashtube. 

We  shall  assume  that  the  radiant  energy  comes  our  uniformly  over 
the  surface  of  the  tube  and  thereby  calculate  values  of  spectral 
radiant  emitta.ice  E0  for  the  lamp.  Thus,  if  the  measured  relative 
energy  output  of  the  lamp  is  G  joules  per  millimicron  per  joule 
input,  and  the  energy  input  is  in  the  form  of  a  nearly  rectangular 
0.75  millisecond  voltage  pulse  containing  f  kilojoules,  the  energy 
output  per  unit  area,  or  spectral  radiant  emittance  E0  is 


joules _  #  f  x  10a  joules  input  ^  1  ,  » 

E°  G  mjx  *  joules  input  (16 .5)ir(l  ,3)cm2  0.75  msec 

The  next  step  toward  obtaining  numerical  values  for  the 
spectral  radiance  N0  is  somewhat  arbitrary;  the  spectral  radiance 
I0  at  the  outermost  surface  of  the  cladding  materials  around  the 
laser  rod  is  taken  numerically  equal  to  the  spectral  radiant 
emittance  E0  of  Eq.  (75)*  This  equality  does  not  imply  that  a 
flashtube  with  a  spectral  radiant  emittance  E0  will  produce,  with 
the  help  of  a  reflecting  enclosure,  such  as  spectral  irradiance, 
but  only  that  it  is  possible  (by  employing  several  i.lashtubes  if 
necessary)  to  produce  an  irradiance  of  the  magnitude  calculated 
from  reference  (27)  by  using  Eq.  (75).  Since  the  radiation  field 
is  to  be  isotropic,  the  connection  between  I0  and  N0  is 
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where  0  is  measured  from  the  normal  to  the  element  of  surface  area. 
It  is  easy  to  verify  that  for  the  case  of  isotropic  radiation. 


TTN0  • 


Thus,  finally, 

No(>0  ,  i Lx  10s.  G  - watts - 

64.35  ir2  mji-cm2  -ste radian 


(77) 


(78) 


The  curves  of  N0/f  which,  under  the  enumerated  assumptions, 
correpond  to  the  spectral  output  curves  for  the  FX-47A  flashtuba  as 
reported  by  Goncz  and  Newell  are  shown  in  Fig.  6.  Values  of  N0  for 
every  5  millimicrons  from  300  to  1000  millimicrons  were  used  in  the 
computer  program. 

4.5  REFRACTIVE  INDICES  AND  ABSORPTION  COEFFICIENTS 

A  neodymium-doped  laser  glass  (#3835  barium  crown  glass) 
containing  4.3  wt$  of  Nd3+  ions  made  by  American  Optical  Company) 
was  selected  for  the  laser  material  of  this  study.  The  required 
refractive  index  values  were  computed  to  five  significant  figures 
from  the  equation 


nia  -  i  +  1.27245  [1  -  0.01/x3  ]-x  -  0.0047655^* 


(79) 


where  is  in  microns.  Measured  values  used  to  ov*ain  Eq.  (79) 
were  the  index  values  1.5259*  1.5219*  1.5198*  1.5-*-fl*  1.5098  and 
1.5093  at  the  wavelengths  0.486  ,  0.5^*  0.589,  O.656,  1.014  and 
1.060  microns  respectively.  The  absorption  coefficient  a^*), 
see  Fig.  7,  was  obtained  from  transmission  measurements  made  on 
plane  parallel  samples  with  a  Cary  Model  14  recording  spectro¬ 
photometer  . 

Tra:  ..  !.ssion  measurements  were  made  with  the  same 
instrument  on  a  samarium  cladding  glass,  EOD-83O,  containing 
9$$  SiHgOg.  The  absorption  coefficient  derived  from  the 
measurements  is  shown  in  Fig.  8a. 

The  condition  of  an  exact  index  match  between  the  samarium 
glass  used  as  a  cladding  and  the  neodymium  glass  of  the  rod  itself 
(to  suppress  whispering  modes)  was  not  duplicated  in  the  computer 
program  because  of  possible  program  halts.  Instead,  the  samarium 
glass  was  given  an  index  na  «  nx  -  0.0001. 
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Figure  6 .  Spectral  radiance  N0  from  the  FX-47A  f lashtube  in  watts 
per  millimicron-cm2 -steradian  divided  by  the  number  f  of  kilojoules 
of  energy  in  an  0.75  millisecond  rectangular  input  voltage  pulse 
to  the  flashtube .  (After  Goncz  and  Newell) 
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Figure  7.  The  absorption  coefficient  a:  in  nun-1  of  AO  5^55 
laser  glass  . 


The  refractive  index  of  the  water  cladding  used  in 
another  of  the  calculations  was  taken  from  published  tables;28  • ii9 
the  absorption  coefficient,  see  Fig.  8b,  was  found  from  trans¬ 
mission  measurements  made  with  a  Cary  Model  14  '■pectrophotometer . 

Pyrex  tubing  was  chosen  as  typical  envelope  material 
for  the  water  cladding.  Refractive  index  values  were  taken  from 
the  relation 

n  »  1.46815  -  0.0080515*2  +  0.0021l34(k2-  0.035)"1 

+  0. 000022788 ( _  0.035)"2,  (80) 

which  was  adjusted  to  fit  measurements  made  on  a  Corning  Pyrex 
#7740  plate .  Transmission  measurements  also  were  taken .  No 
significant  absorption  was  found  except  near  300  millimicrons, 
see  Fig.  8c.  Pyrex  ^740  is  not  an  optical  glass  and  transmission 
characteristics  are  not  controlled,  so  Eq.  (80)  and  the  curve  of 
Fig.  8c  must  be  considered  only  as  representative. 

4 .6  CHECK  CALCULATIONS 

The  first  two  of  the  program  loops  stated  at  the  end  of 
section  4.2  give  an  approximation  to  the  monochromatic  energy 
density  at  a  chosen  point  M  in  the  laser  rod.  The  larger  the 
number  of  rays  traced,  the  better  is  the  approximation  up  to  the 
point  where  machine  round-off  errors  begin  to  interfere.  However, 
the  more  rays  the  machine  has  to  trace,  the  longer  the  computing 
time.  It  is  necessary,  then,  in  starting  the  program,  to  estimate 
the  number  of  rays  to  be  traced  to  the  point  M  for  any  given 
wavelength.  Ideally,  one  would  select  values  for  the  angular 
increments  Act  and  Ary,  compute  an  answer,  then  select  smaller 
values  and  compute  again  —  continuing  the  process  until  the 
answer  remained  constant  to  the  necessary  number  of  significant 
figures . 


An  abbreviated  version  of  this  process  was  carried  out 
for  the  monochromatic  calculations  of  Sooy  and  Stitch30  on  an 
unclad  ruby  rod.  Their  curves  were  duplicated  to  within  ±  3#  by 
taking  Act  =  Ay  ■  5°  .  An  extended  printout  of  the  computation  for 
a  point  near  the  surface  of  the  rod  showed  the  "forbidden  ray” 
phenomenon  explained  in  section  4.3  and  sketched  in  Fig.  5»  The 
minimum  number  of  rays  traced  per  quadrant  was  100,  while  252  rays 
were  traced  when  the  point  M  was  at  the  center  of  the  rod.  The 
greatest  difference  from  the  result  of  Sooy  and  Stitch  occurred 
for  a1R1  »  0.003,  the  least  absorbing  case,  and  x  ■>  0.  The 
calculations  were  repeated  for  that  point  with  Aa  *  Ay  «  2.5°, 
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and  Act  -  5°  ,  Ay  ^  0.625°  which  reduced  the  difference  to  2$  and 
0.5$  respectively.  The  computing  time  for  this  case  and  the  two 
repetitions  were  in  the  ratio  1:4:8.  Estimates  of  the  overall 
computing  time  (a  minimum  of  10  hours  on  an  IBM-360-30  for  ten 
points  within  an  unclad  rod  using  Aa  «  Ay  «  2.5°  for  x  <  0.4r 
and  Aa  =  Ay  -  5°  for  x  &  0.4r,  with  X  =  0.3  (0.005)  1.0  micron) 
made  any  further  refinement  of  the  calculation  inadvisable. 

It  is  believed  that  the  close  agreement  with  the  results 
of  Soov  and  Stitch  verifies  the  computational  method  of  the  first 
two  loops  of  this  program,  particularly  in  view  of  the  considerable 
difference  between  this  method  of  approximation  and  the  method  of 
approximation  used  by  Sooy  and  Stitch. 

The  third  and  fourth  loops  of  the  computer  program  were 
verified  step  by  step  from  an  extended  printout  of  an  actual 
computation  (as  were  the  steps  of  the  first  and  second  loops). 

The  logic  of  the  canputat_  ial  method  was  given  an  overall  check 
by  computing  the  distribution  of  heat  production  in  an  unclad  rod 
and  comparing  the  result  to  similar  calculations  by  Borrelli  and 
Charters ,31  They  calculated  energy  distributions  in  Nd3+  doped 
rods  of  Corning  0580  laser  glass  using  the  spectrum  between  0.3 
and  0.9  microns  obtained  from  an  FT-524  fiashtube  operated  at 
4.4  kv  from  a  400  y,F  condenser  bank.  The  spectrum  was  arbi¬ 
trarily  normalized  against  its  peak  value .  Results  are  presented 
for  several  values  of  the  parameter  pR,  where  p  is  the  concen¬ 
tration  of  doping  ions  in  ions/an3  and  R  is  the  radius  of  the  rod 
in  centimeters.  The  curve  of  rnergy  which  causes  heating  vs.  the 
normalized  radius  interval  when  pR  *  2  x  103°  ions/cm2  has  been 
replotted  in  Fig.  9. 

The  parameters  of  the  check  calculation  verifying  the 
agreement  between  the  method  of  Borrelli  and  Charters  and  the 
method  of  this  report  were  chosen  with  an  eye  on  seme  of  the 
values  used  by  Welling  and  Bickart33  in  their  experimental  work 
on  laser  rods.  Thus,  the  FX-47A  fiashtube  spectrum  corresponding 
to  5300  amp/cm2  was  chosen  as  more  nearly  equivalent  to  the 
output  obtained  by  Borrelli  and  Charters  for  their  FT-524 
fiashtube .  Ifce  AO  3835  laser  glass  has  an  ion  concentration 
of  4.7  x  1020  ions/cm3  so  that  for  a  rod  0.45  cm  in  radius  the 
product  pR  is  2.1  x  1020  ions/cm3. 

The  heating  energy  released  per  unit  volume  in  such  a 
rod  was  calculated  for  11  points  along  the  radius.  For  comparison, 
normalization  at  x  -  0  to  the  value  reported  by  Borrelli  and 
Charters  for  pR  =  2  x  1020  ions/cm2  was  carried  out.  The  agreement 
seen  in  Fig.  9  is  considered  good  evidence  that  the  methods  of  calcu- 
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Normalized  heat  production 


Figure  9.  The  distribution  of  energy  which  causes  heating  in 
unclad  rod  of  laser  glass  with  pR  S  2  x  102  lons/cm  :  ( “H- -  ) 
Corning  0580  glass  (after  Borrelli  and  Charters),  (  O  )  A  5  55 


lation  agree  because  the  physical  properties  of  the  materials 
described  in  each  case  ar-?  quite  similar.  The  slightly  higher 
production  of  heat  energy  in  the  outer  layers  of  the  rod  relative 
to  that  in  the  central  region  as  pP  is  increased  is  predicted  by 
Borrelli  and  Chapters . 

4.7  RESULTS  -  Profiles  of  Heat  Production  and  of  Temperature 

The  spectral  irradiance  at  the  outer  surface  of  the 
laser  rod  and  cladding  assembly,  ns  explained  in  section  4.4,  has 
been  taken  equal  to  the  spectral  radiant  emittance  of  the  FX-47A 
flashtube  as  computed  from  data  reported  by  Goncz  and  Newell. 

These  authors  also  calculated  the  spectral  efficiency  of  the 
lamp  for  the  two  current  densities,  1700  amp/cma  and  5300  amp/cm2, 
and  the  ectral  range,  0.35  to  1.1  microns,  investigated.  By 
taking  the  spectral  efficiency  to  be  the  same  for  the  range  0.3 
to  1.0  microns,  it  is  possible  to  assign  specific  values  not  only 
for  the  spectral  irradiance,  but  also  for  the  irradiance  itself. 

A  current  density  of  1700  amp/cms  is  produced  by  a  square-wave, 

3/4  millisecond  electrical  pulse  of  1  kilojoule.  Assuming  that 
the  fJashtube  is  64.6$  efficient,  646  joules  of  radiant  energy 
appear  in  the  light  pulse .  Thir  is  emitted  uniformly  (by 
assumption)  throuqh  a  surface  area  of  v(l.?)(l6.5)  =  21.45  ir  cm2, 
so  the  radiant  emittance  and  hence  tho  irradiance  corresponding 
io  the  1700  amp/cm2  discharge  is  about  12.8  kilowatts/cm2  over 
the  wavelength  ranya  0.3  to  1.0  microns.  The  irradiance 
corresponding  to  5300  amp/cm2  current  density  is  about  64.1 
kilowatts/cm2  =  (5) (648)/(21 .45  ir)(0.75). 

Four  profiles  were  computed,  all  for  a  pR  value  of 
2.1  x  1C2°  ions/cm2,  or  if  one  prefers,  RF  «  1.9,  where  p  is 
the  concentration  of  doping  ions  in  ions/cm3 ,  F  is  the  ion 
concentration  in  weight  percent,  and  R  is  the  radius  of  the 
laser  rod  in  cm.  The  profiles  of  Fig.  10  shows  the  distribution 
of  heat  production  in  a  Nd3+  doped,  polished,  la^ar  rod  with 
optical  pumping  through  the  side,  (1),  (2)  without  cladding, 

(3)  with  water  cladding  and  (4)  with  samarium  glass  cladding. 

The  curves  have  been  normalized  against  the  irradiance  occurring 
at  the  outer  surface  of  the  l^ser  rod  and  cladding  assembly. 

Curves  (1),  (3)  and  (4)  were  computed  for  the  spectrum  produced 
when  the  current  density  in  che  flashtube  is  1700  amp/cm3 ,  while 
curve  (2)  was  computed  for  a  current  density  of  5300  amp/cm2 . 

Figure  10  has  the  immediate  experimental  interpretation 
that  if  one  were  to  set  up  a  pumping  enclosure  which  could  fill  a 
volume  of  space  with  an  isotiopic  energy  density  of  light  flux 
having  the  rpectral  distribution  of  the  FX-47A  flashtube  or  its 
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Figure  10.  The  distribution  of  heat  production  T  vs  the  normalized 
radius  in  rods  of  AO  3835  laser  glass  with  RF  1.9*  The  ordinate  T 
is  measured  in  watts/cm3  per  watts/cms  irradiance  at  the  surface  of 
che  laser  rod  and  cladding  assembly:  (l)  unclad  rod,  1700  amp/cm2 
lamp  spectrum;  (2)  unclad  rod,  5300  amp/cm3  lamp  spectrum;  (3)  rod 
clad  with  w?.ter  4.17R  thick  (held  in  place  by  an  outer  shell  of 
Pyrex  tubing  0.39R  thick),.  1700  amp/cm2  lamp  spectrum;  (4)  rod 
clad  with  samr.riu"1  glass  0.5R  thick,  1700  amp/cm3  lamp  spectrum 
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equivalent,  then  t'ne  relative  power  dissipation  as  heat  in  a 
laser  rod  with  RF  =  1.9,  (a)  unclad,  (b)  clad  with  water  of 
thickness  4.17  R,  and  (c)  clad  with  samarium  glass  of  thickness 
0.5  R  would  follow  the  curves  (1)  or  (2),  (3)  and  (4)  respectively. 

The  normalization  for  the  distribution  of  heat  pro¬ 
duction  was  adopted  because  it  leads  quiv  y  to  numbers  of  direct 
use  in  laser  pumping.  For  example,  if  the  surface  of  an  unclad 
rod  of  AO  3835  l*ser  glass  with  RF  =  1.9  is  irradiated  by 
1  kilowatt/cm3  o::  light  power  distributed  in  the  spectral  range 
0.3  to  1.0  microns  according  to  the  spectrum  of  the  FX-47A  flash- 
tube  operated  at  1700  amp/cms  current  density,  then  the  rate  of 
release  of  heat  energy  at  the  center  of  the  rod  will  be, 
according  to  curve  (1)  of  Fig.  10,  920  watts/cm3.  By  the  same 
arithmetic,  if  the  light  pulse  from  the  lamp  contained  1  kilo¬ 
joule/cm3,  920  joules/cm3  of  heat  energy  would  be  released. 

It  can  be  inferred  from  curves  (1)  and  (2)  that  the 
change  in  lamp  spectrum  produces  very  little  change  in  the  shape 
of  the  curve,  a  result  which  is  expected  since  the  changes  in  the 
lamp  spectrum  do  not  coincide  with  the  major  absorption  regions 
of  the  laser  glass.  (Compare  Figs.  6  and  7)*  The  higher  concen¬ 
tration  of  heat  energy  in  the  central  region  of  the  unclad  rod, 
due  to  the  focusing  effect,  is  the  predominant  feature  of 
curves  (1)  and  (2l .  When  the  rod  is  clad  with  a  material  of 
index  between  that  of  air  and  that  of  the  rod,  the  focusing 
effect  is  altered  significantly,  curve  (3),  so  that  an  annular 
region  in  the  rod  becomes  the  site  of  the  maximum  heat  production 
per  unit  volume.  Such  a  ring  structure  can  occur  in  an  unclad 
rod  for  higher  concentrations  of  doping  ion  as  shown  by  Borrelli 
and  Charters  and  suggested  by  the  calculations  of  Cooke,  McKenna 
and  Skinner.33  A  hint  that  the  ring  structure  might  occur  in  the 
clad  rod  for  a  lower  concentration  of  doping  ions  can  be  gleaned 
from  the  paper  by  Borrelli  and  Charters,  the  one  by  Cooke  et  al, 
and  one  by  McKenna.34  The  reasoning,  which  is  admittedly  somewhat 
nugatory,  is  as  follows:  Curves  by  Borrelli  and  Charters  indicate 
that  the  distribution  of  heat  production  is  closely  tied  to  the 
distribution  of  absorbed  energy.  Curves  by  Cooke  et  al  indicate 
that  the  distributions  of  absorbed  energy  for  "two-  and  three- 
dimensional"  pumping  are  quite  similar.  Curves  by  McKenna  for 
two-dimensional  pumping  show  that  the  ring  structure  can  be  made 
to  appear  in  a  clad  rod  for  a  much  smaller  concentration  of 
doping  ion  than  is  required  in  the  unclad  rod.  Be  that  as  it 
may,  the  phenomenon  is  undesirable  when  ore  is  trying  to  produce 
a  uniform  distribution.  McKenna  docs  find,  however,  that  the 
thickness  of  the  cladding  affects  the  size  and  shape  of ‘the 
annular  region,  a  fact  which  suggests  that  the  thickness  of 
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the  cladding  should  be  chosen  as  a  function  of  the  doping  ion 
concentration  or,  more  likely,  the  product  pR  in  order  to  produce 
the  most  nearly  uniform  distribution. 

Specifically,  in  Fig.  10,  let  us  compare  curve  (1)  for 
a  polished  unclad  rod  to  curve  (4)  for  a  rod  with  cladding 
thickness  0.5  R.  If  we  nov  imagine  the  clodding  thickness  re¬ 
duced  toward  zero,  we  expecc  the  curve  representing  temperature 
distribution  to  deform  smoothly  from  the  shape  of  curve  (4)  into 
that  of  curve  (1).  At  some  cladding  thickness  the  curve  will 
show  a  minimum  temperature  excursion,  i.e.,  there  should  be  a 
cladding  thickness  which  is  optimum  for  producing  a  uniform 
temperature  distribution. 

The  temperature  rise  in  the  laser  rod  is  nearly  pro¬ 
portional  to  the  heat  production  per  unit  volume .  The  constant 
of  proportionality  for  AO  3835  laser  glass,  with  a  dersity  of 
2.63  gm/cm3  and  a  specific  heat  capacity  Cp  =  0.14  gm-cal/gm-C0 , 
is 

AT(c° )  =  O.65  Em  (joules/cm3).  (81) 

We  turn  again  to  the  FX-47A  flashtube  and  a  rod  of 
AO  3835  laser  glass  such  that  RF  ■  1.9  for  a  >ncrete  example 
of  the  temperature  rise  which  can  occur  in  a  laser  rod.  Suppose 
an  assembly  of  such  tubes  operating  at  a  current  density  of  1700 
amp/cm3  produces  12.8  kilowatts/cm2  irradiance  at  the  surface  of 
the  laser  rod  and  cladding  assembly.  Let  the  duration  of  the 
flash  be  3/4  ms.  Then  9.6  joules/cma  are  incident  on  the  laser 
rod  and  cladding  assembly  per  flash.  According  to  curve  (1)  of 
Fig.  10,  (9-6)(0.92)  =  8.8  joules/cm3  of  heat  will  be  generated 
at  the  center  of  the  unclad  rod  with  corresponding  values  outward 
to  the  edge.  Curve  (1)  of  Fig.  11  shows  the  magnitude  and  the 
distribution  of  temperature  rise  one  flash  of  the  pump  assembly 
will  produce.  Curves  (3)  and  (4)  show  the  temperature  rise  for 
the  same  pumping  when  the  rod.  is  clad  with  water  4.17R  thick  and 
samarium  glass  0.5R  thick  respectively.  Curve  (2)  was  computed 
for  the  unclad  rod  irradiated  by  a  flashlamp  assembly  working  at 
5300  amp/cm£  current  density  and  producing  48.1  joules/cm3  per 
flash  at  the  surface  of  the  laser  rod. 

Experiments  carried  out  by  Welling  and  Bickart  on  rods 
of  AO  3835  laser  glass  provide  a  comparison  between  the  experi¬ 
mentally  observed  distribution  of  temperature  lise  caused  by 
pumping  and  that  predicted  by  the  theory  presented  in  this  report. 
For  RF  as  1.9,  R  =  0.45  cm,  they  report  a  temperature  rise  of  about 
7°C  at  the  center  of  an  unclad  rod  after  one  flash,  wherein  the 
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Figure  11.  Temperature  rise  in  a  laser  rod  of  neodymium  glass 
with  RF  =  1.9  when  pumped  by  a  xenon  flashtube:  (l)  unclad  rod, 
9.6  j/cm2  incident  energy,  (2)  unclad  rod,  48.1  j/cm2  incident 
energy,  (j)  rod  clad  with  water  4.17R  thick,  9.6  j/cm2  incident 
energy,  (4)  rod  clad  with  samarium  glass  0.5R  thick,  9.6  j/cma 
incident  energy. 
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flashlamp  input  energy  to  laser  rod  cross  section  ratio  is 
1000  j/cm3 .  Since  the  7*5  cm  long  rod  had  a  cross  section  of 
O.636  cm3,  the  flashlamp,  probably  of  about  65#  efficiency,  was 
supplied  with  636  joules,  and  radiated  413  joules  per  flash.  If 
the  coupling  between  lamp  and  rod  were  73$  efficient,  301  joules 
reached  the  rod,  or  the  surface  density  of  energy  was  301  /(it) (0.9) 
( T * 5) ■*=  14.2  j/cm3.  By  curve  (1)  of  Fig.  10  and  Eq.  (81),  the 
predicted  temperature  rise  at  the  center  of  the  rod  is  (14.2) 
(0.92) (0.65)  -  8.5cC. 

A  comparison  between  the  experimental  and  the  calculated 
distributions  is  shown  in  Fig.  12.  The  value  of  73$  chosen  for 
the  coupling  efficiency  in  the  preceding  paragraph  was  used  to 
obtain  the  broken  curve  of  Fig.  12  and  represents  a  normalization 
of  the  curve  such  that  the  total  energy  released  as  heat  is  the 
same  for  the  theoretical  distribution  as  it  was  for  the  experi¬ 
mental  one . 

Figure  12  shows  a  noticeable  difference  between  the 
experimental  and  theoretical  distributions  of  temperature  near 
the  surface  of  the  rod.  This  same  phenomenon  has  been  observed 
and  discussed  by  Welling  and  Bickart  for  unclad  ruby  rods .  No 
satisfactory  explanation  of  the  phenomenon  is  known. 
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Figure  12.  Comparison  between  the  temperature  distribution  found 
experimentally  by  Welling  and  Bickart  (solid  curve)  and  the  distri¬ 
bution  calculated  from  the  same  rod  with  similar  pumping  but 
neglecting  internal  reflections  (broken  curve)  . 
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5.  EFFECTS  OF  THE  LASER  FIELD  ON  THE  INDEX  OF  REFRACTION 

OF  THE  LASER  ROD 


5.1  INTRODUCTION 

In  the  following  we  shall  consider  the  effect  of  the 
laser  field  on  the  distribution  across  the  laser  rod  of  changes 
in  index  of  refraction  caused  by  the  field  itself.  In  section  5*2 
we  shall  treat  the  Maxwell  stress  in  the  rod  and  at  the  boundary, 
the  Kerr  effect,  and  the  electros trictive  effect.  In  section  5*5 
we  shall  discuss  the  photoelastic  effect  and  take  as  the  applied 
mechanical  pressure  the  Maxwell  stress  computed  in  section  5.2. 

In  section  5*4  we  shall  attempt  to  justify  the  assumptions  made 
in  the  analysis  and  suggest  a  method  (valid  in  principle)  of 
achieving  greater  laser  beam  uniformity. 

5.2  MAXWELL  STRESS,  KERR  EFFECT  AND  ELECTROS TRICTIVE  EFFECT 
A.  Mechanical  Force 

The  mechanical  force  f  per  unit  volume  due  to  the 
presence  of  an  electric  field  in  a  material  of  dielectric 
constant  e  and  density  p  is  given  by36 

1  *  -  £ E2  9rad  e+  £ 9rad  (E“p  %)  (82) 

It  will  be  our  task  to  find  explicit  expressions  for  e 

and  p  ~  in  the  presence  of  electrostriction  and  any  quadratic 
dp 

effect  in  the  field. 


B.  Electrostriction 


We  wish  to  calculate  the  change  in  density  Ap  owing  to 
the  presence  of  an  electrostrictive  force.  From  the  Helmholtz- 
Lippman  law36  the  change  of  volume  5v  for  an  initial  volume  vQ  is 
given  b  • 


6v 


£v\  1^1  il 

dp  /  8v  3v 


v0E2  9e 
8t  3v 


VqE3  /3e\ 

8*  \3P/E  =  0 


(85) 


T  =  Const. 
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where  p  is  the  pressure 

E  the  electric  field  strength 
e  the  dielectric  constant  of  the  material 
k  the  compressibility. 


The  change  in  density  Ap  and  the  density  p  lead,  using 
Eq.  (83),  to  the  relation 


6v  Ap  E3  / 9  e  \ 

"  v0  “  p  “  Sir  \3P/E  .  0 

(84) 

X  X  | 

I 

T  «  Const . 

So  that  we  have 

n  *- 

*  E3  /  de  \  E3 

ap  -  f  57  il?4 . 0  5  p  5? a- 

(85) 

”  5 

9*  7 

T  *  Const. 

■*  r 

Equation  (85)  expresses  the  electrostrictive  effect  due 

! 

to  the  presence  of  the  field  E. 

- *  ; 

C.  Optical  Kerr  Effect 

-  -  1 

Owing  to  the  presence  of  the  electric  field  E 

the 

**  f 

medium  becomes  birefringent .  The  respective  dielectric 

constants 

?  f 

become37 

i  t  \ 

ee  «  e  +  E3N  (a  -  b)a 

(86) 

-  f 

e0  -  e  +  E3N  (a  +  2b)a 

(87) 

•  ’  j 

where  e  is  the  dielectric  constant  x  to  E 
e 

e0  is  the  dielectric  constant  H  to  E 
and  N,  a,  b,  a  are  constants  of  the  medium. 

Tlie  corresponding  indices  of  refraction  are 


ne  '  V*7  -  'F  [*  +  EatI  (a'b)  1?]* 

(88) 

n„  -  ^  -  yjT  [1  +  E»N  (a+2b)  £•]«  n, 

(89) 

50 


Mtli'-i- 


In  the  Kerr  effect  measurements,  both  DC  and  optical, 
one  measures  the  quantity 

no  “  ne  *  E?N  2^  -  biXe2  (90) 

where  Bl  is  the  Kerr  constant. 

In  our  case  and  for  our  purpose  we  are  interested  in 
the  absolute  change  of  index 

n0  -  n  =  nE2N  (a  +  2b)  —  (91) 

It  will  be  noticed  that  in  the  measurement  of  Bx  only 
the  coefficient  b  in  Eq.  (90)  is  determined.  In  the  case  of 
Eq.  (91)  the  coefficient  a  has  to  be  determined  independently. 
Using  Eq.  (87)  we  can  write 

e0  -  e  «  E2N  (a  +  2b)a  e  B  §1  (92) 

8tt 

where  we  have  defined  another  constant  B  which  we  shall  subse¬ 
quently  use . 


D.  Conditions  at  the  Boundary 

(1)  Surface  Force  per  Unit  Area 

In  this  section  we  wish  to  calculate  the  Maxwell 
stresses  at  the  boundary  of  the  laser  rod.  We  shall  have  to 
transform  the  force  per  unit  volume  of  Eq.  (82)  into  a  force 
per  unit  area.  The  first  term  in  Eq.  (82)  becomes  (omitting  the 
cumbersome  factor  of  l/8ir  for  the  moment) 

-  E3  grad  e  =  grad  (-  e  E3)  -  e  grad  (-  E2)  (93) 

The  second  term  on  the  right-hand  side  will  mostly 
concern  us  since  it  is  not  readily  amenable  to  the  desired  trans¬ 
formation.  From  vector  analysis36  we  have 

grad  (a* b)  =  (A*v)  B  +  (B-y)  A  +  A  X  (v  X  b)  +  B  X  (v  X  a)  (94) 

We  shall  henceforth  only  handle  the  electrostatic  case 
(see  justification  in  section  5-4),  hen^e 

V  X  E  =  0. 

Then  we  have, 

e  grad  E2  =  2e  (e*v)  E  =  2(d*v)  E  (95) 
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where  D  «  eE  is  the  electric  displacement  vector.  In  order  to 
make  headway  with  the  term  D*v  in  general  (as  opposed  to  cartesian 
coordinates)  we  shall  have  to  draw  more  heavily  on  vector  analysis. 
We  define  the  indefinite  or  dyadic  product  of  two  vectors  A  and  B 
as  A;B.  Taking  a  third  vector  C  we  define  the  scalar  product  of 
the  dyad  A; B  by 

(A;B)  •  C  -  A(B*C)  (96) 

and 

c  •  (A;B)  -  b(c-a)  (97) 

Thus  the  product  i  vector,  and  is  different  depending  on 
whether  C  follows  >  recedes  the  dyad.  For  the  operator  nabla  v 
we  define  for  a  gi  ector  F 

v;F  «  lim  -  /  n ■  F  dS  (98) 

V-0V*/(S) 

where  n  is  the  normal  vector  to  a  closed  surface  S  enclosing  a 
volume  V.  This  definition  is  analogous  to  the  definition  cf  the 
gradient  of  a  function  0,  which  looks  as  follows 

7 0  «  lim  —  f  n  <t>  dS  =  grad  0  (99) 

M'-'ls) 

Similarly  for  the  divergence  of  an  operator  we  have  the 

definition 

v  •  F  =  lim  I  f  n  •  F  dS  =  div  F  (100) 

V-.0  V^(S)  x 

We  can  now  find  a  clearer  meaning  to  the  operator  A  •  v 
as  given  in  Eq.  (9^). 

Let  A  be  a  vector  which  is  not  to  be  varied  in  the 
limiting  process  given  by  Eq.  (98),  then 

lim  I  /  (A  •  n)  FdS  -  A  *  lim  —  f  (n;F)dS  »  A*(v;F)  (101) 
V-OV^(s)  |_ V-*0  V*'  (cj)  J 

and  from  Eq.  (97)  get  for  the  last  term  in  Eq.  (101) 

A* (7  j  F)  =  (A*  V )  F  (102) 
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Given  two  vectors  a  and  b  we  wish  to  calculate  the 
divergence  of  the  dyad  a;b.  We  have 

v* (a,b)  =v*(a;b0)  +v‘(a0;b)  (103) 

where  the  suffix  zero  indicates  that  the  corresponding  quantity 
is  not  to  be  varied  when  applying  the  operator  v. 

Applying  Eq.  (100)  to  the  first  term  on  the  right  of 
Eq.  (103)  yields 

V- (a;b0)  =  iim  !  f  n‘(a;b0)dS  (104) 

v-0  V  ‘'(s) 

Using  Eq.  (97)  this  is  equal  to 

limi  f  (n-a)b.  dS  =  lim  ^  [  (n-a)  dS  b„  -  (va)b  (105) 

v-o  v-^  (s )  L  vJ(s)  J 

This  follows  from  Eq.  (100)  and  we  have  dropped  the  suffix  zero 
on  the  extreme  right  hand  side. 

So  that  we  have 

v*  (a;b0)  -  (v-a)b  (io6) 

Let  us  do  the  same  thing  for  the  second  term  on  the 
right  of  Eq.  (103).  Using  Eq.  (100)  and  Eq .  (97)  we  have 

V-  (an b)  «*  lim  I  f  n*  (a0  ;b)  «  lim  i  f  n*(aD)bdS  (107) 
v  01  '  v.ovi(s)  V-0  vJ  (g) 

This  last  expression  on  the  right  leads  by  Eq .  (101) 
and  Eq .  (102)  to 

lim  i  f  (n*a0)b  dS  «=  (a*v)b  (108) 

V-°  vj(s) 

So  that  we  have 

v*(a0;b)  *  (a*v)b  (109) 

Assembling  Eqs.  (109)  and  (106)  we  get  the  final  result 
for  Eq.  (103) 

v-  (a;b)  =  (v*  a)b  +  (a*v)b  (H°) 
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We  now  derive  what  may  be  called  the  tensor  form  of 
Gauss'  Theorem. 


(V) 


(a;b)  dr  -  f  n. (a;b) 
J(S) 


dS 


(111) 


written 


L 


(v) 


written 


Using  Eq.  (110)  the  left  hand  side  of  Eq.  (Ill)  can  be 

r-  (a*b)  dt  ■  f  (v*a)b  dr  +  f  (a*v)b  dT  (112) 

•/(V)  J{v) 

Using  Eq.  (97)  the  right  hand  side  of  Eq.  (Ill)  can  be 


C  n-(a;b)  dS  -  f  (n*a)b 
7(S)  J 


dS 


(113) 


So  that  the  tensor  form  of  Gauss'  Theorem  becomes,  using 
Eqs.  (112)  and  (113) 

y(v*a)b  dr  +  f  (a*v)b  dr  =  f  (n«a)b  dS 

(v)  */(v)  J(t) 

By  invoking  Eq.  ( 114)  we  can  now  write 
f  |"e(v«d)  +  (d*v)e|  dr  -  f  E(n-D)dS 

J{V)  L  J  J(s) 


(114) 


(115) 


In  our  case,  since  we  do  not  have  any  free  charges  in 
the  laser  rod  we  have  v*D  ■  0,  hence  Eq.  (115)  becomes. 


f  (d* V )E  dT  *»  f  E(n*D) 
•/(v)  J(r\ 


dS 


(s) 


(1x6) 


Equation  (116)  provides  us  with  the  sought  for  relation  enabling 
us  to  pass  from  a  volume  integral  of  Eq.  (95)  to  a  surface  integral. 
Integrating  Eq.  (99)  over  a  volume  V  we  have 

f  v  0  dT  -  f  grad  $  dx  ■  f  n  $  dS  (117) 
*'(v)  •'  ( s\ 


f(v) 


(S) 


f  We  are  now  in  a  position  to  transform  the  volume  integral 

}  (v)f  dr  of  Eq.  (82)  into  a  surface  integral.  Using  Eq.  (117) , 
the  volume  integral  of  the  second  term  of  Eq.  (82)  becomes 

h  /vl9rad(ESp  dF)dT-5F  /slnESf>  3FdS  (118) 
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Similarly,  the  volume  integral  of  the  first  term  on  the 
right  hand  side  of  Eq.  (93)  becomes 

'  U1*  <eEM  dC  ’  '  5?  /  ,n  dS  (119) 


The  volume  integral  of  Eq.  (95)  ba comes 

&  f(v)  E  9rad  E°  dT  -  Wr  f.g)  E(n-D’  “  <12°) 

where  use  has  been  made  of  Eq.  (116). 

Assembling  Eqs.  (118),  (119)  and  (120),  the  volume 
integral  of  Eq.  (82)  becomes 

/  £  at  -  -L  /  E(n-D)  as  -X  f  nEa  It  -  o  (121'i 

J(v)  It  “'(s)  8t  J(s)  (  p  ip)  1 ’ 


So  that  we  have  a  surface  force  t 


t  -  ~  D(E*n)  -  IT- (*  ~  P 


(122) 


This  is  the  Maxwell  stress . 


(2)  Maxwell  Stress  at  the  Eoundary 

Since  we  wish  to  calculate  the  Maxwell  stress  at  the 
boundary,  we  shall  use  Eq.  (121)  and  integrate  over  a  surface 
given  by  Fig.  13.  The  boundary  separating  medium  1^  and  2.  is 
given  by  BB'.  We  assume  an  area  dS  perpendicular  to  the  plane 
of  the  figure  and  indicated  by  a'0'  for  medium  2  and  our  area  dS 
indicated  by  afi  for  medium  1_.  The  area  which  shows  in  the  plane 
of  the  figure  as  a'a  and  0'|3  will  be  chosen  so  small  that  its 
contribution  to  the  surface  integral  will  be  neglected. 

Hence,  we  have 

f  f  dx  ■  dS 

Av) 

-  & [E’2  (E  -  ^  p)a  'V(e'tp)])  (12J) 


1 5J-  [e,Ea(Ea.n)  -  e1B1(Es.n)j 
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We  shall  assume  that  medium  2  is  the  vacuum.  Hence,  the  pressure 
p  at  the  boundary  is  given  by 


p  -  57  [E.  (®.  •  »)  -  <Ei  '  ")]  -  £  ( 

Ea2 -Ej  2 

oj  (124) 

where  we  have  put  (  e  -  —  pi  =  a. 

\  dp  A 

(123) 

Let  us  set  up  a  coordinate  system  where  the  field  E  is 
decomposed  into  wo  components i  En  parallel  to  the  vector  n,  and 
Ep  parallel  to  tv.<  boundary.  Then  we  have 

Ea  *  n  =  Ean  5  Ei  *  n  *  Em 

(126a) 

V  =  Eln  +  Elp  ;  Ei3  =  Ein  +  E?p 

(126b) 

E  =  1  E  +  1  E 
n  n  p  p 

(126c) 

where  ln  and  lp  are  unit  vectors  respectively  parallel  and  perpen¬ 
dicular  to  the  vector  normal  n.  Inserting  Eqs.  (126a),  (126b)  and 
(126c)  into  Eq.  (125)  yields. 


P‘w[(1n 


Ean  +  Ip  Eap 


) 


E 


an 


E 


in 


+  1 


P 


(127) 


The  boundary  values  of  the  field  are 

ei  Ein  =  Ean  (128a) 

Eap  -  Elp  (128b) 


Inserting  Eqs.  (128a)  and  (128b)  into  Eq.  (127)  yields, 
P  =  W  [(Xn  El  Ein  +  1pEip)  El  Ein  “  (Xn  El«  +  *p  ^p)  Einj 
-  £  Ein  +  Eip)  +  °  ("‘?n  +  E?p)  (129) 
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We  note  that  the  terras 

*p  [ei  EipEin  "  ei  EipEin]  in  E<3 •  (129)  cancel. 

We  are  thus  left  with  normal  stresses  only.  Rewriting  Eq.  (129) 
yields 

p-T^-fl  e,aE?_-e,l  E?  1 
4ir  L  n  1  ln  1n1nJ 

-  &  («.a  +  E?p)  +  t  “  (E?n  +  E?p) 

or  finally 

P  a  brn  [E'n  ^  "  2£l  +  a*  +  E?p  (a  "  *)]  (W) 

Deperding  on  the  sign  of  p,  we  shall  have  a  surface 
pressure  or  tension  at  the  boundary. 

E.  Calculation  of  f 

We  shall  derive  ein  expression  for  f  as  given  by  Eq.  (82) 
taxing  electrostriction  and  Kerr  effect  into  account.  Using 
Eq.  (85)  which  yields  the  expression  for  the  change  in  density 
due  to  electrostriction ,  the  total  density  p  is  given  by 

p  »  p0  +  Ap  -  p0  +  p0  a  «  p0  (l  +  aj  (131) 

Using  Eqs.  (92)  and  (85) ,  the  total  dielectric  constant  resulting 
from  electrostriction  and  Kerr  effect  becomes, 

e  -  e0  +  Ac  -  Ap  +  B  (132) 

It  should  be  noted  that  B  in  Eq.  (132)  can  also  be  made 
to  include  any  quadratic  contribution  to  the  dielectric  constant 
other  than  the  Kerr  effect. 

It  will  be  our  task  now  to  use  Eqs.  (131)  and  (132)  to 

find  a  more  explicit  form  for  Eq.  (82).  We  start  with  the  second 

term  on  the  right  hand  side  of  Eq.  (82)  first. 


58 


Putting 


(133a) 

(133b) 


=  p  =  x  +  xy 


(133c) 


We  have 


Hence 


Thus 


3e  _  3e_  9ja  ^ 

3x  3u  3x  3u  +  ^ 


3e_  =  1  .  3e 

3u  1  +  y  3x 


u  3e  _  u  3e 
3u  1  +  y  3x 


/  x  +  xy\  3£  a  x  de 
[l+yldx  3x 


(134) 


(135) 


(136) 


Reverting  back  to  the  physical  values  as  given  by 
Eqs .  (133a),  (133b),  (133c),  Eg.  (136)  becomes. 


3  £  3e 

P  7~  ■  Po  an 

dp  3p0 


(137) 


„  de  -  v 
Po  3Po  ' 


(138) 


is  a  quantity  which  can  be  determined  from  photoelastic  measurements 
as  will  be  seen  below.  In  our  case  y  is  the  value  of  p  3je  before 

dp 

the  application  of  the  electric  field.  In  this  work  we  shall 
henceforth  assume  that  y  is  constant  in  the  laser  material  and 
does  not  vary  with  position  in  the  medium. 

Before  turning  to  the  first  term  on  the  right  hand  side 
of  Eq.  (82),  we  wish  to  rewrite  Eq.  (132)  after  substituting 
Eq.  (131)  for  Ap.  Then  Eq.  (132)  becomes. 


e  =  e0  + 


p0  a+  B  *  +y|ia+Bp- 

8ir  8u  Sir  8ir 


=  30  +  |1  (ya  +  B) 


(139) 
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ji  . . .  nil' 


And  after  putting 


b  s-  y  a  +  B 

Eq.  (139)  can  be  written  as 


p2 

e°  +  b  S7 


(140) 


(141) 


Then  the  first  term  on  the  right  hand  side  of  Eq.  (82) 
becomes,  assuming  e  uniform  throughout  the  laser  rod,  and  using 
Eq.  (l4l) 


-  ^  Ea  grad  e  E2  gr?id||i  (ya  +  b) 


B 


b  grad 


(142) 


The  second  term  or.  the  right  hand  side  of  Eq.  (82) 
becomes,  using  Eq.  (137) 


Or  3rad 


(E3(%)  ’  orart  (Y  1?) 


Nov/,  since  we  have  assumed  y  uniform  throughout  the 
laser  rod  we  can  Lake  out  y  from  under  the  grad  operation  on  the 
right  hand  side.  Bence,  Eq.  (143)  yields 


grad  (e*  p  ^ 


"  V  9rad  IS 


(144) 


Assembling  Eqs.  (142)  and  (144),  we  get  the  final  form 
for  Eq.  (82),  namely 


-  b  grad 


+  y  grad 


E2 

8ir 


(145) 


F.  Calculation  of  y  in  the  Clausius-Mossotti  Case 


ii  i 

**  1 


*  3 

I  I 


i  1 

i  :  i 


Since  for  a  large  class  of  glasses  it  is  reasonable  to 
assume  that  the  Clausius-Mossotti  relation  applies,  it  will  prove 

useful  to  calculate  Hp  for  this  case .  The  Clausius- 

Y  =  P7 

Mossotti  relation*  dP 


*See,  for  example,  M.  Born,  Optik  (Springer,  Berlin,  1933)  P*  339 
or  Born  and  Wolf,  Principles  of  Optics  (Pergamon  Press,  1959)  P.  86 
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. . .  .  . * . 


"•tfllflJVtWtttttVM 


can  be  written  for  our  purposes  as 


1 


e  -  1 


-  c  P 


e  +  2 

where  C  is  a  constant 
By  ordinary 


de 

dp 


and  p  is  the  density, 
differentiation 

-  l][c  +  2] 


(146) 


(147) 


Hence,  we  get  the  desired  relation 

pff=3“C=  -  !][  =  +  21  s  V  (118) 

The  change  Ae  in  e  due  to  electrostriction  or  changes 
in  density  Ap  becomes 

(fle)p  ”  ^  *P  -  J  0  -  1]  0  +  2]  f  (1*9) 

The  change  Ae  in  e  due  to  the  Kerr  effect  or  any  other 
change  in  polarizability  Aa  becomes 

'  s>a  -§**  *  -  !]  [*  +23  f  (150) 

When  the  Clausius-Mossotti  law  applies,  we  get  the 
interesting  relation 


P 


de 

dp 


°  \  [e  -  1]  [e  +  2] 


(151) 


Here  a  is  buried  in  the  constant  C  which  itself  becomes  a  variable. 
The  constant  c  is  proportional  to  the  polarizability  o. 

Since  the  index  of  refraction  n  is  related  to  the 
dielectric  constant  e  by  the  relation  e  =  n3  we  can  also 
express  y  in  terms  of  n  as 


or 


—  =  2n  p 
dp  P 


P 


dn  =  7 
dp  ~  2n 


dn 

dp 


(152) 

(155) 
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G .  Maxwell  Stress  at  the  Boundary 


(^)  In  the  case  for  which  the  Clausius-Mossotti 
relation  holds. 

We  shall  evaluate  the  coefficients  Efn  and  E2 
in  Eq.  (1J0).  Using  Eq.  (151)  we  have  " 

-  2 ex  +  a  ■  exa  -  2ea  +  e1  -  L  (2e1a  +  2e^  -  4) 

6 

-  g-  [^eia  "  8ei  +  4  ]  -  g-  [ei3  -  2ei  +  1]  “  3  [*»  "  <154) 

Similarly  for  the  coefficient  of  E2p  we  have 

a  -  1  -  -  1  j^8  +  -  4  J  -1  -  £-2e,8  +  4e,  -  2  j  i 

-  -  |  [V  -  26,  +  1  ]  -  -  |  (e,  -  l)8  (155) 

Assembling  both  terms  Eq.  (130)  becomes 


p  ‘5F  (e'  '  1)3  3  [2  E?n  '  e'p]  ' 


(156) 


We  notice  that  the  first  term  is  a  tension  force  and 
the  socond  term  a  pressure  force . 

As  an  illustration,  we  shall  calculate  the  pressure  at 
the  edqe  of  the  cylinder  in  the  case  of  a  normal  component  Ein 
only  being  present.  From  Eq.  (156)  we  then  have 

P  -  |  O  -  lf^-Efn  (157) 

which  is  a  tension. 

For  1  Mw  per  cm2  power  we  have 

2 

3  (1.52  -  l)2  •  106  •  107  _  .  ,  „  _ 

p  =  - - - - - =  350  dyn/cro2  a*  0 .36  gm/cm2  (157a) 

>  *  1™ 
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. . .  . * . •■••• . .  .IHtliljj; 


since 


Ef_  ^  W 
8ir  c 


where  W  =  power/cm2 

c  =  velocity  of  light. 

For  200  Mw/cm2  power  we  have 

p  .»  7  *  101  dyn/cm2 


(15Tb) 


This  indicates  that  there  is  a  mechanical  force  being  applied  on 
the  cylinder  as  v  result  of  the  presence  of  the  E  field.  Tfriis 
force  will  contribute  to  a  change  in  index  dee  to  the  photoelastic 
effect  as  we  ^hall  see  in  section  5«3« 

(b)  In  the  general  case . 

We  look  first  at  the  coefficients  5  =  e,2  -  2ex  +  a 
and  a  -  1  of  Eg.  (130).  Using  Eqs.  (139)  and  (140)  we  have  for  ex 


£i  ■  *  -  £o  +  fp  * 


a«eD  -  Y+  fi-b 


(158) 


(15S0 


For  ease  of  notation  let  us  put 


b  =  Ae 


(160) 


Then  we  shall  lead  up  to  the  evaluation  of  6 
ex2  -  (e0  +  Ae)2  =  e03  +  2Ae  e0 


(161) 


where  we  have  neglected  the  term  (Ae)2 . 


We  now  have  for  5 


e12  -  2ex  +  a  =  ecs  -  ec  +  Ae  (2e0  -  1)  -y 


(162) 


If  we  assume  the  largest  value  for  y  to  be  3*^0  as 
computed  in  Table  II  for  the  Pocke Is’  glasses  and  putting  n  *  1 .5# 
we  get  for  Eq.  (162) 


5  =  -  0.60  +  3.5^b. 


(163) 
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II 


Values  of  Y  for  Various  Glasses 


S205 

0428 

0658 

02154 

01571 

0500 

S57 

Observed  p  ~ 
3p 

0.358 

(0.106) 

0.392 

0.483 

0.512 

0.600 

0 .865 

n 

1.508 

1.51c 

1.545 

1  -570 

1.644 

1.751 

1.96 

2n 

3.016 

3.024 

3.090 

3.14 

3.288 

3.502 

3.92 

*  ■  f>  fr 

1.08 

1.21 

1.55 

1.62 

2.1 

3.40 

op 


This  corresponds  to  a  change  in  index  of  refraction 
An  *=  ,  so  that  An  due  to  y  a  ~  is,  (n  ■  1.5) 

An  -  £§.  «  ULL  IP"6.  -  6  •  lO-7  (170) 

2n  3.0 

For  the  constant  B  in  Eg.  ( 1 40 ) ,  ve  shaxl  assume  that 
it  is  of  the  same  order  of  magnitude  as  the  Kerr  constant.  In 
Tauern39  we  find  for  the  absolute  Kerr  constant  in  glass >  the 
values 

1.7  •  10-13 
5.9  *  10"13 
5.8  •  10-13 
8.1  •  lO”13 

The  coefficient  a  -  1,  becomes 

q  ..  J  B  g  _  y  _  1  +  ^*—  b 

0  8 T 

Using  y  *  1 .21  as  given  in  Table  II 
a  -  1  *  0.04  +  ^  b 

We  therefore  note  that  only  if 

eoa  ”  ec0  +  7  ) 

>  is  closr  to  zero 

e0  -  Y  “  1  J 

the  term  b  can  be  important . 

or 

However,  let  us  assume  B  ■  0(1),  and  a  laser  power  of 
200  Mw/cm2 .  Then 

P  =  n  6  -Ss.  =  200  ■'  -  6.6  •  104  dynes/cm3  (175) 

or  5  •  1010 

This  pressure  at  the  boundary  will  produce  later  a  change  in  index 
due  to  the  photoelastic  effect. 


(171) 

(172) 
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The  value  of  Eq.  (173)  should  be  compared  to  that  of 
Eq.  (157b).  It  is  of  the  same  order  of  magnitude. 

We  terminate  this  section  with  Table  II  where  we  have 
calculated  the  values  of  y  for  several  glasses  as  given  by  Pockels40 
from  his  observed  values  of  gn 

P  a?  ' 

We  have  simply  used  Eq.  (153)  and  from  some  inferences  in  the 
text,  we  have  found  the  corresponding  indices  of  refraction  which 
ve  are  taking  as  given  by  Pockels. 


5.3  PHOTOELASTIC  EFFECT 

In  this  section,  we  wish  to  apply  the  results  of  our 
previous  work.  From  the  Maxwell  stresses  calculated  in  the 
previous  sections  it  is  now  possible  to  evaluate  their  effect  on 
the  index  of  refraction.  We  shall  not  redo  here  the  theory  of  the 
photoelastic  effect.  Instead  we  shall  refer  the  reader  to  Adams 
and  Williamson41  (A  &  W) .  We  shall  use  their  notation.  We  use 
A  &  WT  s  Fig.  3,  which  becomes  our  Fig.  1*4.  The  mechanical  stress 
is  applied  along  OY.  We  are  interested  in  the  absolute  changes 
of  index  ny  -  n  for  light  polarized  in  the  y  direction,  and  in 
nz  -  n  for  light  polarized  in  the  z  d.rection.  We  reproduce  in 
Table  III,  A  &  W s  Tab-e  3.  We  reproduce  here  the  formulas  of 
interest  for  this  section.  They  are 


JL - -  *  JL  (  %.  -  (A  &  W  Eq.  5) 

n  2R  \  v  v  / 

JL -  .  -  ( -2  5^  +  ^)  (A  &  W  Eq.  3B) 

r.  E  \  v  v  / 

Y 

which  we  approximate  to 

„  -  n  -  -  (-8.E4) 
y  e  \  v  v/ 

and  finally 

"  I  [<!  -  0)  v  -  0  v]  (A  4  ”  Eq>  3C) 


(174) 


(175) 
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Figure  14.  Notation  used  by  Adams  and  Williamson  —  drawing  to 
accompany  elementary  discussion  of  optical  effects  of  stress. 

The  thrust  P  is  applied  in  the  direction  OY.  The  ray  of  light 
enters  the  cube  of  glass  in  the  direction  OX,  becomes  elliptically 
polarized  and  is  treated  as  two  rays  vibrating,  respectively,  in 
the  directions  OY  and  OZ.  Ordinarily  the  ray  vibrating  along  OY 
travels  with  the  higher  velocity,  that  is,  ordinary  glass  under 
unidirectional  compression  behaves  like  an  optically  negative 
uniaxial  crystal . 
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TABLE  III 

Results  of  Pockels'  Measurements  on  Optical  Effects  of  Stress* 
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which  we  approximate  to 


P 

n  -  n  =  n  — 
z  E 


-  a)  2.  _  a 
v 


(176) 


The  birefringence  quoted  xa  Table  III  corresponds  to 
Eq.  (174).  We  have  computed  :*or  the  same  glasses  (identified  by 
their  Pockels  number)  the  values  as  given  by  Eqs.  (175)  and  (176) 
respectively,  for  P  *  1  kg/cm2  .  This  is  indicated  in  Table  IV. 


TABLE  IV 

Index  Changes  for  1  kg/cm2  Pressure 


Pockels 

n  (°) 

Sample  No. 

y 

205 

5.06 

428 

-8,38 

658 

1.06 

2154 

1.99 

1571 

3.44 

500 

4.76 

57 

9.45 

n 

n 

z 

n 

10”  8 

48.30 

•  10“8 

10- 8 

19.30 

•  10-8 

10- 7 

4.83 

•  10" 7 

10- 7 

4.89 

•  10-7 

10-7 

6.08 

•  10- 7 

10- 7 

6.13 

•  10" 7 

10- 7 

7.52 

•  10-7 

It  is  interesting  to  note  that  although  the  composition  of  all 
glasses  listed  (except  sample  57  which  is  positively  birefringent) 
are  negatively  birefringent,  glass  sample  428  exhibits  a  negative 
value  for  ny  -  n.  For  boundary  pressures  of  6.6  •  104  dynes/cm3 
as  given  by  Eq.  (173)  the  change  in  index  An  will  be 


6.6  •  10* 

9.80  •  106 


or 

nt°)  -  n 
z 


(1  kg  weight  =  9.80  •  106  dynes) 


(177) 


which  would  correspond  for  glass  sample  57  and  n  -  n 
to  Y 


An 


6.6  •  104  -  9.45  •  IQ"7 

9.80  •  106 


6.38 


10-8 


9.45  *  10-7 
(378) 
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which  is  small.  Tie  question  of  the  variation  of  An  across  the 
laser  rod  can  be  dealt  with  adequately  from  the  foregoing  work, 
once  the  field  distribution  across  the  rod  is  given. 


Vfe  wish  to  terminate  this  section  by  pointing  out  that 
it  is  possible  to  express  £n  (and  thus  y)  in  terms  of  the 

P  3p 

photoelastic  constants  p  and  q,  as  follows. 


P 


Bn 

dp 


This  was  first  derived  by  Pockels.40 


(179) 


5.4  DISCUSSION 

In  all  the  foregoing  we  have  treated  the  electrostatic 
case  only.  The  physical  justification  resides  in  the  fact  that  p 
cannot  follow  the  field  at  optical  frequencies.  The  DC  component, 
however,  will  be  present  and  exercise  the  influence  we  have 
described.  If  the  description  is  adequate  it  would  follow  that 
it  ought  to  be  possible  in  principle  to  compensate  index  changes 
across  the  laser  rod  (leading  either  to  convergent  or  divergent 
beams)  by  applying  an  electrostatic  field  with  the  proper  field 
distribution  (which  can  also  be  made  to  vary  in  time ) .  We  have 
also  assumed  equivalence  between  optical  and  DC  Kerr  effect, 
which  is  reasonable  for  gless. 


5.5  SUMMARY 

In  sections  5.2  to  5.5  we  have  treated  the  problem  of 
the  influence  of  the  laser  field  itself  on  the  index  of  refraction 
of  the  laser  rod  (with  special  interest  in  a  glass  rod) .  We  have 
developed  the  inter-relationship  between  Maxwell  stresses  (both  in 
the  bulk  material  and  at  the  boundary),  electrostriction,  photo¬ 
elastic  effect  and  Kerr  effect  (both  optical  and  DC  in  the  case 
of  glass).  We  have  developed  the  formalism  to  a  point  where  the 
distribution  across  the  laser  rod  of  changes  in  index  of  refraction 
can  be  readily  computed  provided  one  is  given  the  distribution  of 
the  laser  field  in  the  rod. 
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6.  END  REGION  STRESS  ANALYSIS 


The  presence  of  a  non-uniform  temperature  distributior 
in  a  circular  laser  cavity  gives  rise  to  two  significant  effects 
which  will  modify  the  optical  path  of  light  passing  through. 

First,  the  strain  patterns  which  accompany  a  temperature  gradient 
serve  to  alter  the  optical  pathlength  for  different  radial  positions 
within  the  rod.  Where  the  material  shows  a  stress-optic  behavior, 
the  index  for  the  two  polarizations  of  light  will  also  be  different. 
Second,  as  there  is  also  a  finite  curvature  of  the  end  surface 
which  accompanies  a  non-un-form  distribution,  there  will  be  a 
discrete  bending  or  refraction  of  the  ray  as  it  passes  through 
the  surface . 

In  the  present  effort,  a  program  was  initiated  to 
evaluate  the  stresses  and  the  optical  distortion  they  produce 
in  a  finite  length  cavity  under  axially  symmetrical  temperature 
distribution.  Both  temperature  and  stress  fields  thus  are 
assumed  to  exhibit  no  dependence  on  the  azimuthal  coordinate,  6. 

In  practice,  these  conditions  can  be  closely  approximated  with 
isotropic  pumping.  For  a  laser  rod,  then,  in  which  pumping 
intensity  is  uniform  along  its  length,  and  temperature  variation 
within  the  rod  is  an  arbitrary  function  of  radius,  the  stress 
field  may  be  resolved  exactly  over  the  major  portion  using  plar.a 
strain  theory.  The  region  in  which  the  plane  strain  formulation 
will  not  be  valid  is  at  the  extreme  ends,  where  according  to 
St.  Venant’ s  principle,  self-equilibrating  mechanisms  will  be 
set  up  in  order  that  specific  boundary  requirements  are  met. 

Since  measurable  changes  occur  in  the  stress  field  within  one 
two  diameters  of  the  rod  end,  giving  rise  as  well  to  surface 
rmation,  cne  would  expect  related  changes  in  optical  path 
differing  from  a  plane  strain  analysis  based  on  total  length. 

A  detailed  study  of  the  end  region  stress  field,  therefore,  was 
undertaken . 

Treatments  of  the  end  problem  in  cylinders  raported  in 
the  literature  are  meager  because  the  solution  long  withstood 
attempts  even  for  an  approximate  development.  In  principle, 
however,  a  solution  can  be  formulated  exactly  by  relating 
stresses  and  displacements  to  the  Love  displacement  function, 

L(r,z)  as  described  in  Ref.  42.  Now,  if  one  set?  the  condition 
that  L(r,z)  be  biharmonic,  i.e.,  satisfies  the  equation 

v4 (l)  =  v8 vs (l)  =  0  ( 180) 
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where 


=  £L_  +  i  i.  +  a 1_ 

3r3  r  3r  9za 


r  and  z  being  cylindrical  coordinates ,  the  equations  of 
equilibrium  and  compatibility  for  axially  symmetric  systems  will 
then  be  satisfied.  Further,  if  the  boundary  conditions  to  be 
satisfied  are  no  forces  exerted  on  the  curved  surface,  and  self- 
equilibrating  symmetric  normal  and  shear  tractions  on  the  end 
corresponding  to 


a  m  o,  t  ■  0 
r  * 


at  r  a  a 


J  az(r)rdr  =0  at  z  =»  0 


the  (l)  function  will  take  the  form 


L(r,z)  -  e^Vfyr)  + 


2(l-v)  +  YJ0(Y)/Ji  (y) 


e”^ZYr  ji  {yr )  (181) 


where  here  y  becomes  a  solution  to  the  eigenvalue  equation 


J%(y)  m  +  2(l-v) 

j?(y)  y 2 


(182) 


In  the  absence  of  external  forces  but  where  loading  is  caused  by 
thermal  expansions,  the  derivations  hold  if  the  temperature  function 
can  be  reduced  to  a  symmetrical  system  of  forces  distributed  at 
the  end. 


While  the  above  equations  do  constitute  a  system  of 
functions  for  the  exact  solution  of  the  problem,  a  difficulty 
in  use  arises  from  the  fact  that  the  roots  of  Eq.  (182)  consist 
of  real  and  imaginary  parts .  The  task  of  calculating  these 
complex  roots  has  been  found  to  be  exceedingly  great,  and  the 
real  and  imaginary  parts  of  Eq.  (181)  give  rise  to  increased 
mathematical  difficulty  at  the  boundary.  Therefore,  the 
expressions  described  do  not  provide  practical  workable 
solutions . 


In  an  attempt  to  circumvent  some  of  these  difficulties, 
Horvay43  has  introduced  an  approximate  variational  method  of 
solution  in  which  two  Sadowsky-Sternberg44  stress  functions, 
given  a  product  representation,  0(r,z)  =  f(r)g(z),  and 
4>(r,z)  =  F(r)G(z),  are  utilized  to  determine  the  stress 
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distribution.  Briefly,  by  introducing  the  stress  functions  in 
terms  of  stresses  into  the  expression  for  complementary  strain 
energy  and  taking  the  variation,  6U,  equal  to  zero  there  will 
result  Euler  equations  with  constant  coefficients  from  which  the 
axial  functions  g(z)  and  g(z)  can  be  determined.  g(z)  and  G(z) 
will  be  established  as  exponentially  decaying  sinusoidal  functions 
of  tht  form, 


“GZ 

e  (cos0z  +  £L  sin0z)  . 

3 

The  a’ s  and  0' s  are  real  and  imaginary  components  to  the  complex 
root  y  =  a  +  i0  of  the  eigenvalue  equation 

Ay4  -  By2  +  C  =  0 

where  A,  B,  and  C  are  constants  taken  from  the  Euler  equations 
resulting  above . 

The  functions  f(r)  ana  F(r)  are  radial  polynomials 
which  must  satisfy  certain  boundary  conditions  and  rule  of 
orthogonalization .  They  must  be  selected  so  that  with  their 
appropriate  derivatives,  they  will  constitute  a  system  of 
boundary  tractions,  ar,  ae,  02  and  Trz  at  z  =  0,  which  will 
closely  approximate  the  set  of  end  tractions  imposed  by  the 
problem.  The  variational  method  permits  the  end  tractions  to 
be  expanded  employing  conventional  Fourier  expansion  techniques 
which  alleviates  much  of  the  mathematical  difficulty  incurred 
at  the  boundary  with  the  exact  method.  Once  successful  in 
determining  matching  polynomials  and  their  derivatives,  one 
may  evaluate  the  Euler  constants  A,  B,  and  C  from  which  axial 
functions  g(z)  and  G(z)  will  follow.  To  calculate  the  product 
stress  functions  and  resultant  principal  stresses  which  they 
represent  then  becomes  straightforward  using  available  Sadowsky- 
Sternberg  relationships.  While  primary  obstacles  inherent  with 
an  exact  formulation  can  be  removed  by  this  approach,  it  does 
remain  that  considerable  mathematical  involvement  and  judgment 
is  still  required  for  execution  of  the  method. 

Additional  papers  (though  more  distantly  related  to 
the  present  problem)  dealing  with  analytical  solutions  to  the 
axial  symmetrical  loading  of  the  cylinders  of  finite  length  have 
appeared.  In  Refs.  45  and  46,  stress  analyses  have  been  con¬ 
ducted  on  the  isothermal  cylinder  where  prescribed  tractions  and 
displacements  on  the  lateral  surface  together  with  the  ends  was 
investigated.  The  problem  of  thermal  stresses  in  ho  ll-cv  cylinders 
with  fairly  small  length  to  radius  ratios  is  taken  up  by  Bellamy.47 


75 


The  method  of  analysis  here  is  based  on  a  complementary  strain 
energy  theorem  for  thermoe.lasticity  and  approximate  solutions  for 
a  cylinder  of  particular  proportions  have  been  worked  out. 

In  our  present  work,  we  have  undertaken  to  handle  the 
thermal  stress  problem,  initially  at  least,  by  way  of  a  numerical 
method  of  solution.  There  appeared  to  be  these  special  advantages 
in  this  approach  over  the  analytical  method  outlined  above: 

(1)  The  governing  equations  for  a  numerical  relaxation 
stress  analysis  could  be  developed  with  minimum  difficulty  and 
application  to  the  end  region  of  a  cylinder  seemed  straightforward. 

(2)  The  equations  could  be  programmed  for  solution  on 
a  high  speed  computer  allowing  investigation  of  many  different 
temperature  distributions,  including  axial  temperature  variation, 
without  constituting  a  new  boundary  value  problem. 

For  such  axially  symmetrical  systems,  we  have  equations 

by  Hoyle40 

I3x  ■  0 

£!*  „  ia*  .  SSL.  r  —  (183) 

0z3  1-v  dr 

which  have  been  found  suited  to  relaxation  solutions.  The  stress 
variables  r,z)  and  x(r,z)  take  the  form 

j  .  r 

dr 


x 

1-v  3r 

where  0  and  Q  are  arbitrary  functions  of  the  coordinates,  r  and 
z.  E,  a,  v  and  T(r,z)  above  are  Young’s  modulus,  thermal  expansion 
coefficients,  Poisson’ s  ratio  and  temperature  at  a  specific  co¬ 
ordinate  location,  respectively.  Is  is  the  differential  operator: 

If _ I 

dr3  r  dr  dz3 
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The  principal  stresses,  ar,  ae,  etc.  can  be  expressed  in  terms 
of  the  functions  and  x  as 


^  +  (l-v)x  1  3  , .  ,  »  EaT 

a  m  L - — L - — -  +  x)  -  - — 

r  r2  r  3r  1-v 


o  -  i  ii  .  t  -  -  i  2*. 
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(184a) 
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where  it  may  be  verified  by  substitution  that  the  equilibrium 
equations  given  by 
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are  identically  satisfied. 


The  bounds  of  the  system  impose  certain  other  requirements 
on  the  equations  as  well.  Since  the  surfaces  of  a  laser  rod  will 
remain  essentially  free  of  external  force,  the  boundary  conditions 
to  be  applied  take  the  form 
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(186) 


where  here  *1*  and  TmT  are  direction  cosines  of  the  normal,  N,  at 
boundary  points  under  examination .  It  may  be  observed  that  on  the 
horizontal  surface,  these  expressions  reduce  to 
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Figure  15-  Point  designations  for  relaxation  grid. 


Figure  15  represents  a  schematic  representation  of  the 
axial  cross-sectional  portion  of  a  laser  rod  used  for  analysis. 

The  cylinder  axis  becomes  a  'lound  of  symmetry;  also,  at  the  left, 
about  1.5  diameters  from  the  end,  another  arbitrary  bound  of 
symmetry  is  established  along  a  section  where  stress  patterns 
no  longer  exhibit  dependence  on  the  axial  coordinate,  z.  The 
end  region  so  defined  between  these  bounds  and  the  free  surfaces 
is  then  subdivided  into  a  convenient  number  of  zones  or  points  to 
form  the  relaxation  grid.  A  general  point  in  this  array  is 
assigned  the  designation  ij,  referring  to  the  row  and  column 
number.  In  conventional  notation,  a  neighbor  to  the  left 
becomes;  i,  j-1;  one  to  the  right;  i,  j+1;  and  so  on. 

Appropri  he  dependence  between  points  in  the  chosen 
array  is  established  by  expressing  the  governing  expressions  in 
finite  difference  form.  Using  the  calculus  finite  differences. 
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first  and  second  order  differentials  may  be  expanded  by  the 
following  good  approximations 


where  a  is  the  grid  spacing . 


(187) 


Working  equations  can  now  be  set  up  embodying  the  stress 
functions,  and  x,  vrhich  will  be  representative  of  each  category 
of  points  in  the  system.  Since  we  are  dealing  with  two  functions, 
we  have  in  general  two  equations  to  be  served  simultaneously  at 
array  points.  By  introducing  Eqs.  (187)  into  Eqs.  (I85)  1  iveming 
the  interior  region  there  results  the  following  two  expressions 
which  are  suitable  for  numerical  analysis 


+  -  W 


(189) 


Points  which  are  located  along  the  horizontal  curved  boundary 
must  obey  formulas  ( 188)  and  (189)  and  satisfy  boundary  conditions 
(186)  as  well.  Also,  tl*2  fact  that  positions  i+l,j  lying  above 
the  boundary  plane  are  fictitious  requires  suitable  substitution 
for  the  variables  at  these  points.  By  solving  Eq.  (188)  for 
x  i+i  i  and  (189)  for  ^i+l,j  and  writing  a  third  equation  from 
(184a)  at  the  cylindrical  surface,  p. 
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we  may  combine  Eqs .  (188),  (189)  and  (190)  to  eliminate  xi+1 
and  fi+l,j  from  the  boundary  equation.  In  addition,  since 
Ti+l,j  is  also  a  fictitious  temperature,  the  temperature 
difference  (Ti+i^j-  Ti_i#j)  appearing  in  Eq.  (190)  must  be 
expressed  in  terms  of  real  temperatures.  For  this  purpose,  a 
MacLaurin  series  expansion  is  used  to  expand  the  temperature 
differential  about  surface  point,  ij,  in  the  following  manner 


If  we  now  express  these  terms  as  finite  differences,  we  find  the 
above  reduces  to 

(“)„■  k  (3Tii +  Ti-2,j  - 

which  is  equivalent  to 

*Ti+l,j  "  Ti-1 , j ^ 


By  making  this  substitucion  and  performing  the  indicated  algebra 
the  boundary  equation  becomes  in  final  form 
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For  points  along  the  vertical  boundary  we  find  from  con¬ 
dition  (186)  that  the  axial  stress  goes  to  zero  and  may  therefore 
write  from  (184b) 


a 
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i+1  >  j 
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(192) 


This  requires  that  ^i+l,j  equal  ^i-l,j  and  implies  that  be  a 
constant  on  the  boundary.  Hence,  by  combining  Eqs .  (188),  (189) 
and  (192)  to  eliminate  Xj.^j+:L  and  V'i-j+l*  the  final  version  of 
the  boundary  equation  appears  as 


*ij  + 
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(193) 


For  points  lying  on  an  interior  bound  of  symmetry, 
additional  conditions  will  hold  and  the  applicable  equations 
above  would  be  modified  where  possible.  Thereby,  at  the  left 
radial  bound. 


Xi,j-1  =  Xi,j+1*»  ^i,j-l  *  ^i,j+l 


For  the  cylinder  axis,  where  i  =  1,  it  results  that  equations 
(188)  and  (189)  become  indeterminate  because  of  a  division  by 
zero  in  the  second  and  fourth  terms.  Therefore,  at  axial 
boundary  positions  a  somewhat  specialized  mathematical  formulation 
is  required.  Beginning  with  the  first  of  the  equilibrium  Eqs.(l85), 
we  find  for  r  =  0  this  reduces  to 
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O'  T 


rz 


0 


on  the  axis. 


79 


I 


3y  combining  with  Eqs .  (l84a,d)  and  performing  the  necessary 
algebra,  the  above  expression  becomes 

where,  for  the  equation  to  be  satisfied  when  r  -  0, 

[*i.3  +  (1-v)  *1.1]  -  °- 

From  the  system  of  interdependent  equations  now  defined, 
the  values  of  x^j  and  for  each  point  is,  in  principle,  determined 
knowing  specifically  only  the  temperature  distribution  imposed  on 
the  network.  The  solution,  however,  demands  that  the  equations 
become  simultaneously  satisfied  which,  in  effect,  requires  that 
residual  quantities  Q^j  and  R^j  be  reduced  to  identically  zero  in 
all  equations . 

Solving  this  problem  by  hand  relaxation  methods  would 
be  prohibitive,  even  for  a  small  array,  because  of  the  tedious  and 
lengthy  computation  involved.  Therefore,  a  considerable  portion 
of  the  task  was  devoted  to  preparing  a  computer  program  which 
would  perform  the  relaxation  operation  automatically. 

A  brief  description  of  the  reduction  technique  developed 
follows.  Handling  up  to  a  300  point  array,  the  program  computes 
and  stores  the  magnitude  of  the  residual  quantities,  Q^j,  and  R^j, 
for  all  points  based  on  selected  values  of  Xj,  ^  and  which  are 
initially  guessed.  The  relaxation  procedure  involves  scanning  the 
stored  arrays  for  the  highest  value  (positive  or  negative)  of  Q 
or  R  which  exists  and  modifying  the  x  and  ^  values  for  the 
corresponding  point  by  a  computed  increment  according  to 

1 

1 
i 

(194) 


where  x{-j  and  become  new  values  ~->r  the  stress  function 
associated  with  the  point,  ij.  The  relaxation  control  variables, 
X  and  0,  may  take  on  values  between  0  and  2  as  prescribed  by  the 
degree  of  under-relaxation  or  over -relaxation  required.  Upon 
altering  values  of  Xij  and  ,  it  obviously  necessitates  re¬ 
computing  the  magnitude  of  residuals  at  neighboring  points  which 
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are  affected.  By  this  process  of  selection  and  correction,  the 
iterative  procedure  is  continued  until  all  residuals  are  diminished 
to  within  a  specified  range  from  zero. 

Using  then  the  stress  function  values  which  have  been 
found,  the  principal  stresses  are  computed  employing  formulas  (184). 
The  corresponding  thermally  induced  axial  strain  field  is  obtained 
using  a  three-dimensional  stress  strain  relationship  of  the  form 

e  =  —  [a  -  v(o  +  a  )1  +  aT  (^95) 

z  e  L  z  r  e  J 

And  the  axial  displacement  of  the  end  face  due  to  variation  of  the 
strain  field  over  the  end  region  is  given  by 

w(r)  =  f  5P  ez(r)dz 

*'0 

where  the  integration  is  performed  over  1.5  diameters. 


The  change  in  optical  pathlength  for  a  ray  of  laser 
light  traversing  the  cavity  end  region  may  now  be  calculated 
using  the  pertinent  expressions  developed  by  Quelle,49  and 
modifying  them  slightly  so  as  to  include  the  axial  variation 
of  strain.  The  resulting  expressions  become  for  a  ray  polarized 

radially 
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and  for  one  polarized  tangentially 
APe(r)  =  (n-1) ez(r)dz 
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where  n  is  th*i  optical  index  of  the  glass;  p  and  q,  optical 
constants  to  be  determined  experimentally;  V0  ,  the  velocity  of 
light  in  the  unstressed  medium;  and  Bx  and  By,  the  characteristic 
stress  optic  coefficients.  To  obtain  the  total  change  in  pathlength 
for  the  optical  cavity,  requires  that  the  above  results  be  combined 
with  additional  results  for  the  portion  of  the  rod  in  which  plane 
strain  applies. 

The  computer  provides  the  principal  stresses  and  strains 
for  each  position  in  the  relaxation  array  as  superposed  on  the  end 
portion  of  the  red.  To  determine  the  axial  optical  path  change, 
then,  through  this  region  requires  that  a  numerical  point-by-point 
integration  be  made  of  the  separate  quantities  defined  in  Eqs.  (196) 
and  (197).  It  was  found  expedient  to  perform  this  integration  by 
hand  using  a  desk  calculator. 

In  the  initial  testing  of  the  computer  program,  con¬ 
siderable  difficulty  was  encountered  in  getting  selected  problems 
to  converge .  The  principal  source  of  difficulty  was  subsequently 
determined  to  be  the  result  of  choosing  too  few  radial  subdivisions 
in  the  array  in  an  attempt  to  minimize  computer  operating  time 
during  the  evaluation  runs.  When  at  least  eight  radial  divisions 
were  designated  for  an  arbitrarily  selected  2.54  cm  diameter  rod, 
the  residual  quantities,  Qij  and  Rj.j>  would  tend  to  zero  and  the 
problem  converged  properly.  It  occurred,  similarly,  that  a 
higher  limit  also  existed  for  the  number  of  radial  points  chosen 
for  which  convergence  could  be  assured.  For  the  work  done  the 
number,  nr,  of  radial  grid  divisions  fell  within  the  limits 

n 

16  —  <  40 

a 

where  a  is  the  radius  of  the  rod  in  inches.  For  a  number  of  cases 
which  fell  outside  these  limits  convergence  was  not  obtained. 

The  relaxation  coefficients,  X  and  <f>,  were  found  tc  play 
governing  roles  in  the  speed  to  which  relaxation  could  be  achieved. 
Ihe  most  satisfactory  value  for  both  coefficients  in  this  work 
seemed  to  be  unity  inasmuch  as  under -relaxation  tended  to  make 
computing  too  slow,  and  over -relaxation  caused  the  solution  to 
diverge  at.  certain  stages  due  presumably  to  inherent  limitations 
on  the  mathematical  process . 

To  verify  the  accuracy  of  the  computer  solutions,  it  was 
sufficient  to  show  that  the  results  obtained  compared  favorably 
with  known  solutions  in  two  limiting  cases.  Specifically,  for  the 
long  rod,  we  kr^»  that  the  stress  distribution  through  a  cross 
section  beyond  >out  1.5  diameters  from  the  end  must  be  given  by 
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a  plane  strain  solution.  In  addition,  for  the  infinitely  thin 
element,  the  disc,  the  results  obtained  must  agree  with  the  plane 
stress  solution  which  is  common  to  all  radial  sections  .  Figures  16 
and  IT  provide  an  illustration  of  how  the  computer  results  compart 
with  hand-calculated  results  using  plane  strain  and  plane  stress 
theory  in  these  respective  configurations.  First,  a  2.^4  cm 
diameter  AO  3835  neodymium  glass  rod,  was  considered  with  a  11.1°C 
linear  temperature  rise  from  center  to  edge .  Figure  16  plots  the 
principal  stresses  machine  computed  for  a  cross  section  1 .€  diameter 
from  the  end  using  a  minimum  network  array  of  9  x  27  points. 

Figure  17,  similarly,  plots  the  computed  stress  pattern  at  the 
midplane  of  a  2.40  cm  diameter  glass  disc  using  a  21  x  3  array 
and  a  2.22°C  linear  temperature  rise  between  center  and  edge. 
Diameter  to  thickness  ratio  for  the  disc  was  10  to  1 .  Negligible 
change  was  indicated  in  the  stress  distribution  across  the  disc 
and  axial  stress  remained  essentially  zero.  Treating,  in  addition, 
the  case  of  a  short  cylinder  with  a  diameter  to  length  ratio  of  2, 
it  was  found  that  an  intermediate  solution  to  plane  strain  and 
plane  stress  occurred  at  the  midplane . 

Fairly  good  agreement,  therefore,  exists  between  the 
computed  results  and  the  exact  theory  for  these  cases.  The  gradual 
deviation  indicated  in  the  principal  stresses  with  increasing 
proximity  to  the  axis  is  due  to  small  residual  errors  in  the 
computed  stress  functions  which  magnify  in  turn  as  r3  tends  to 
zero  in  Eqs.  (l84a,d).  Since  a  finite  error  is  introduced  from 
the  transition  of  differential  Eq.  (I83)  to  finite  difference 
form,  small  discrepancies  between  results  cannot  be  removed 
completely.  Error  te;.*ms  generated  in  using  a  square  mesh  array 
are  generally  of  the  order  of  a2  where  f  a*  is  the  spacing  of 
the  points.60  Thus,  the  magnitude  of  these  terms  can  be  reduced 
by  going  to  more  divisions  and  thereby  decreasing  the  mesh  size  . 

The  benefit  of  this  is  demonstrated  in  Fig.  17  where  21  radial 
points  on  the  disc  were  chosen  and  considerably  improved  agreement 
was  achieved  near  the  axis . 

Stress -strain  patterns  in  the  end  region  which  contribute 
to  an  axial  optical  path  variation  are  of  interest.  Accordingly, 
the  simple  arbitrary  case  described  above  using  a  11.1°C  linear 
temperature  gradient  in  a  2.54  cm  rod  is  treated  further.  Computed 
results  showing  the  change  in  principal  stresses  and  axial  strain 
near  the  end  are  depicted  in  Figs.  18  and  19  for  radial  positions 
chosen  at  .477  and  .954  cm,  respectively,  from  the  axis.  Clearly, 
the  solution  agrees  with  the  plane  strain  results  beyond  about 
1  diameter  from  the  end  but  changes  significantly  in  the  end  region 
as  predicted.  Contrary  to  previous  reasoning,  though,  the  stress 
distribution  at  the  end  surface  of  the  rod  does  not  approach  the 
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Figure  16 .  Comparison  of  computer  solution  and  plane  strain  theory 
for  a  cross  section  about  1.6  diameters  from  end  of  rod. 


STRESS  (kg/cm2 ) 

Compression 


Figure  1?.  Comparison  of  computer  solution  and  plane  stress 
theory  for  midsection  of  disc. 
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Figure  18.  stress-strain  distribution  at  0.477  cm.  Radial 
distance  from  axis  in  a  2.54  cm  diameter  0855  laser  rod  with 
11.1°C  temperature  difference  between  center  and  edge. 
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Figure  19.  Stress-strain  distribution  at  0.95^  cm.  Radial 
distance  from  axis  in  a  2.5^  cm  diameter  08j5  laser  rod  with 
11.1°C  temperature  difference  between  center  and  edge. 


ra^  il  id  tangential  results  for  the  disc,  which  from  theory 
would  v  .pproximate ly  4/5  of  the  plane  strain  results  shown. 

Using  computed  stress  and  strain  results,  the  axial 
pathlength  at  various  radial  positions  was  calculated  using 
Simpson' s  rule  and  integrating  Eq.  (196)  above,  ior  radially 
polarized  light,  over  the  end  region.  The  equation  accounts  for 
three  effects;  namely,  the  stress  optical  change,  the  index  change 
with  temperature,  and  the  elastic  deformation  or  strain.  The 
integrated  result  cf  the  latter  component,  the  axial  strain,  gives 
the  physical  displacement  resulting  in  the  end  surface.  Illustrating 
their  relative  contributions  to  the  total  path  change  in  AO  3&35 
glass,  these  components  are  plotted  separately  in  Fig.  20  along 
with  tleir  summated  result.  The  path  difference  or  distortion  is 
expressed  relative  to  the  axial  position.  It  is  interesting  to 
note  that  for  a  linear  temperature  gradient,  a  nearly  linear 
relationship  is  preserved  in  all  components.  This  suggests  the 
possibility  of  compensating  for  differential  path  changes  at  the 
ends  in  such  a  specialized  case  by  utilizing  a  glass  having  an 
appropriate  thermal  index  coefficient  value,  an,  for  tie  rod 
length  involved  where  from  above 


For  comparison,  the  optical  path  change  in  the  end  region  was 
computed  also  for  this  temperature  distribution  usir  lane  strain 
theory  and  the  results  were  higher  by  about  7#  f-om  che  computer 
solution  presented.  In  a  1  meter  long  rod,  it  turns  out  this 
error  would  amount  to  about  1 

A  similar  optical  path  change  can  be  determined  for 
tangentially  polarized  light  by  evaluating  the  various  components 
in  accordance  with  Eq.  (197)  above.  For  AO  5^55  laser  glass,  the 
optical  coefficients  B|(  and  Bx,  were  measured  to  be  6.1  and  8.2 
Brewsters,  respectively.  The  optical  path  for  the  radial 
polarization  is  calculated  to  be  about  13 #  greater  than  that  for 
the  tangential  polarization  at  the  point  of  maximum  deviation 
which  occurs  at  the  surface  of  the  rod.  Since  radial  and 
tangential  stresses  become  equal  at  the  center  of  the  rod,  it 
follows  that  differences  here  due  to  polarization  diminish  to 
zero. 


In  this  study,  attempt.®  were  made  to  assess  end  effects 
for  two  other  radial  temperature  patterns  which  could  be  considered 
more  typical  of  conditions  set  up  during  l-.^ing.  The  first 
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Figure  20.  Optical  path  variation  for  the  end  region  of  a  2.5^cro 
AO  3835  neodymium  glass  rod  with  a  linear  center  to  edge  temperature 
gradient  showing  relative  influence  of  the  thermal  index  change, 
axial  displacement  and  stress  retardation  effect. 


temperature  curve  was  derived  from  experimental  gain  distribution 
measurements,  obtained  under  another  program  in  using  a  2.54  cm 
neodymium  AO  3855  roughened  surface  rod  pumped  at  23.5  kj  input. 
The  profile  of  the  second  curve  was  established  from  the  computer 
calculation  performed  on  this  contract  for  a  polished  0.9  cm 
diameter  rod  of  the  same  glass  {Fig.  11,  curve  2).  Both 
temperature  patterns  represent  significant  non-linearity 
with  radius . 

In  generating  the  first  curve,  temperature  rise  was 
considf  ced  to  be  proportional  to  the  increase  in  gain  coefficient, 
0,  across  the  rod.  This  relationship  proves  basically  valid  since 
the  inversion  of  ions  in  a  four  level  laser  material  can  be  shown 
to  be  relatively  linear  with  the  absorption  of  incident  photons 
per  volume  element.  The  specific  equations  become 

0  = 

L 

AT  cc  3 

where  G  is  the  measured  gain;  L  the  length  of  the  rod,  and  AT, 
the  temperature  rise  of  the  corresponding  position  for  which  0 
is  calculated.  The  proportionality  between  AT  and  0  must  be 
determined  for  the  pumping  geometry  considered. 

The  fraction  of  electrical  energy  absorbed  in  the  laser 
cavity  as  heat  was  actually  investigated  experimentally  by  making 
calorimetric  measurements  on  12.9  mm  and  18  mm  rods  of  AO  3835 
glass.  Each  rod  was  fired  at  an  input  level  of  9500  joules  usin<7 
two  linear  flashlamps  close -coupled  to  the  rod  in  a  conventional 
reflector  configuration  without  cooling.  Immediately  following 
the  pulse,  the  rod  was  removed  from  the  system  and  placed  into  a 
cylindrical  Cewar  flask  where  heat  could  be  exchanged  with  a 
prescribed  quantity  of  water  contained  within  the  flask.  Heat 
production  in  the  rod  was  then  determined  by  making  careful 
temperature  measurements  on  the  water  initially  and  after 
equilibrium  conditions  were  reached.  A  Leeds  and  Northrup 
precision  potentiometer  was  employed  to  record  temperature. 

The  results  are  described  in  Fig.  21.  In  one  case,  the 
flashlamps  were  jacketed  with  pyrex  tubing  which  served  to  absorb 
selectively  much  of  the  incident  UV  and  IR  radiation;  in  the  other 
case  the  flashtubes  were  non-jacketed,  and  this  accounts  for  the 
greater  heat  production.  One  point  is  included  also  from  earlier 
work  in  which  the  heat  generated  in  a  2.54  cm  rod  was  determined 
on  a  fully  assembled  system  where  cooling  water  was  circulated 
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Figure  21.  Experimental  results  for  heat  production  in  AO  5^55 
glass  laser  rods  pumped  in  a  cross-coupled  silver  reflector  wrap. 
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over  the  rod  within  a  pyrsx  jacket.  In  this  investigation, 
differentia),  water  temperature  measurements  were  recorded  in 
time  between  the  inlet  and  outlet  positions  to  the  cavity.  The 
area  represented  by  curve,  then,  was  numerically  integrated  to 
establish  the  heat  quantity  released  from  the  rod.  Fairly  good 
continuity  of  the  data  points  was  found.  Therefore,  it  could  be 
estimated  with  some  certainty  that  between  25  and  30$  of  the 
energy  appeared  as  heat  in  the  2.54  cm  water  clad  rod  under  study. 
Using  the  30$  level  and  the  gain  information  supplied  on  this  rod, 
a  temperature  diffe  nee  of  about  2.8°C  between  center  and  edge 
is  calculated. 

Computed  results  for  the  optical  path  variation  in  the 
end  legion  and  surface  contour  change  are  plotted  in  Fig.  22  along 
with  the  temperature  profile  derived.  The  optical  path,  as 
illustrated,  was  averaged  for  the  radial  and  tangential  polari¬ 
zations  for  purposes  of  comparing  this  curve  with  the  athermalized 
case  where  temper ature  effects  would  be  compensated  for  by  proper 
regulation  of  an.  As  suggested  in  section  3 •  to  obtain 
athermalization  in  hire fr ingen c  glasses  in  a  practical  way,  the 
light  must  oass  first  through  the  laser  in  one  polarization  and 
return  through  the  other.  an,  then,  must  be  averaged  from  the 
stress  optical  effects  in  the  two  orientations.  The  theory 
presumes,  however,  that  the  stress  distribution  will  be  axially 
uniform  between  the  cavity  bounds.  In  a  rod  of  finite  length, 
the  stress  distribution  changes  considerably  from  plane  strain 
near  the  ends,  which  may  warrant  that  conditions  for  athermalization 
be  re-examined. 

The  significance  of  the  ends  is  apparent  from  Fig.  22. 
Actual  path  changes  over  this  region  obtained  by  the  computer 
method  is  represented  by  curve  C  and  the  results  from  another 
calculation  invoking  plane  strain  analysis  throughout  is  shown 
by  curve  D.  Since  the  expression  for  curve  A,  is  a  member  of 
the  composite  expression  for  both  curves  C  and  D,  changes  can  be 
implemented  in  the  latter  curves  simply  through  control  of  the 
temperature  coefficient  of  index  values,  0^,  in  A.  Theie  is  a 
proportionality  evidenced  here  between  curves  D  and  A  and  this 
makes  it  possible  to  compensate  for  the  temperature  effect  com¬ 
pletely  in  a  long  rod  by  making  sufficiently  negative  as 
described  earlier .  Curve  D  would  then  coincide  with  the  vertical 
axis  of  the  graph.  The  difference  between  curves  C  anc'  D,  however, 
would  remain  essentially  fixed;  this  being  the  finite  optical  path 
distortion  remaining  in  an  " athermalized"  glass  rod  due  to  the  ends. 
For  the  present  temperature  distribution,  the  path  difference 
between  center  and  edge  is  about  0.3  microns  (for  two  ends)  when 
athermalized.  This  results  in  very  slight  negative  power  for  the 
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overall  rod  under  t»  se  athermalized  conditions.  The  variation  in 
surface  contour  is  described  by  curve  E. 

To  correct  for  end  effects,  the  value  of  should  be 
chosen  to  be  less  negative,  so  that  curve  C  rather  than  curve  D 
lies  along  the  axis.  For  the  length  shown  (i.e.,  1.25  diameters) 
the  corrected  value  should  have  the  approximate  value 
Oj{  =  .72  a^.  Here  0^  is  the  value  which  will  bring  curve  D  to 
the  axis.  The  value  of  is  that  which  will  bring  curve  C  to 
the  axis  at  rp  «  .75. 

However,  for  a  longer  length  rod  the  end  effect  becomes 
relatively  less  important.  At  a  length  of  1  m,  the  corrected 
value  a^,  will  be 


2.5^ 


50  -1.25  ‘ 
50 


2.54 
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.98  a 


A  slightly  more  complex  temperature  distribution  is 
treated  by  the  second  example .  Figure  23  depicts  the  results 
of  the  computer  solution  for  this  case  and  shows  s^:...swhat  more 
perturbation  due  to  the  ends.  Plane  strain  results  and  actual 
optical  path  change  were  plotted  here  for  the  average  of  radial 
and  tangential  polarizations .  It  may  be  observed  that  as  much 
as  1-1/4  microns  path  difference  exists  between  the  actual  and 
the  plane  strain  curves  at  about  60#  of  the  radial  distance  to 
the  edge.  As  with  the  previous  case,  this  amounts  to  residual 
distortion  that  would  be  present  having  employed  an  athermalized 
glass.  For  this  case,  however,  the  end  face  displacement, 
represented  by  curve  B,  shows  considerably  greater  irregularity. 


In  conclusion,  it  was  demonstrated  that  for  laser 
glasses  which  are  not  designed  to  compensate  for  thermal  stress 
effects,  the  differences  in  optical  pathlength  near  the  ends 
could  become  virtually  obscured  when  considering  the  cumulative 
pathlength  variation  generated  within  the  long  rod.  The  ends 
would  account  for  a  few  percent  change,  at  most,  from  the 
analytical  result  based  cn  plane  strain  calculation  for  the 
full  length.  On  the  other  hand,  though,  it  was  shown  for 
athermalized  glasses  that  end  effects  in  the  rod  could  be  the 
principal  source  of  distortion  left.  Since  path  changes  ir  the 
end  region  are  not,  in  general,  related  directly  to  the  temperature 
distribution,  compensation  for  these  effects  would  require  special 
modifications  beyond  regulation  of  the  glass  properties  alone. 
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Figure  23.  Relative  temperature  change,  end  face  distortion,  and 
optical  path  variation  for  the  end  regiCn  of  an  0.9  cm  diaiti.  3835 
glass  laser  rod  pumped  at  153  j/cm2  . 
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The  most  suitable  scheme  to  use  would  depend  to  seme 
degree  on  geometry  em*  d  and  temperature  pattern  imposed. 
Theoretically,  only  the  thin  disc  configuration  has  axially 
uniform  stress.  Effective  athermalization  of  a  series  of  these 
elements  would  then  be  acnieved  by  utilizing  a  glass  with  the 
appropriate  value  and  correcting  for  the  axial  surface 
displacement  in  use  by  immersion  in  a  suitable  index  ir  itching 
fluid.  For  the  rod  configuration,  where  a  significant 
temperature  structure  occurs,  the  addition  of  properly  designed 
correction  optics  to  the  output  system  may  substantially  reduce 
the  end  problem. 
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